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ABSTRACT 



Large raindrops (greater than 2.2 nun diameter) that strike a water surface at 
terminal velocity are capable of creating bubbles that radiate significant underwater 
acoustical energy. Previous studies have revealed a positive correlation between 
underwater sound spectral levels during rainfall and the number of large raindrops 
present. Therefore, laboratory' measurements have been made of the underwater 
sound generated by large raindrops. Using the laboratory' measurements, smoothed 
energy density spectra for various sizes of large raindrops are determined. These 
spectra are then used to compute a predicted underwater sound spectrum due to 
rainfall for rainfall rates of 15 mm/hr and 100 mm/hr, assuming an exponential 
(Marshall-Palmer) raindrop size distribution. The resulting spectra are compared to 
underwater sound spectra measured at sea during periods with similar rainfall rates. 
The predicted rainfall spectra are comparable to the measured rainfall spectra. 
Possible reasons for differences are discussed. An inversion technique for obtaining 
the raindrop size distribution from the rainfall acoustical spectrum is presented. An 
alternate approach for obtaining the required inversion matrix is suggested for future 



work. 
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I. INTRODUCTION 



Precipitation is an important factor in air/ocean climate modeling and weather 
forecasting. However, measurement of precipitation in the open ocean is a difficult 
venture. Because of platform motion and disturbance of the airflow pattern in the 
vicinity of the platform, rain gauge measurements at sea are unreliable. Moreover, 
rain gauges measure rainfall onl> at the host platform, and so spatial variability in the 
precipitation cannot be determined. Indirect methods forrainfall measurement include 
radar backscatter attenuation measurements and analysis of satellite imagery'. These 
methods tend to be rather inaccurate. An independent means of verifying the rainfall 
measurements obtained by remote sen.sing is therefore desirable. 

The underwater sound produced by rainfall striking the ocean surface may 
provide another means of measuring rainfall rates and raindrop size di.stributions at 
sea. The shape of the underwater sound spectra in the presence of rain is much 
different from the spectral .shape when rainfall is absent (Lemon and Farmer, 1984). 
Also, the .spectral levels when rain is present can be many decibels higher than the 
spectral levels due to wind generated noise, as .shown in Figure 1. Thus, the under- 
water noi.se generated by rainfall is identifiable, which suggests that a correlation may 
exist between the precipitation rate and underw-ater spectral levels. 
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Figure 1. Spectral characteristics of ocean noise (Nystuen and Farmer, 1989). 

However, attempts to develop direct empirical relationships between rainfall rate 
and underwater sound intensities have only been moderately successful. First, sig- 
nificant surface wind effects on the rainfall underwater acoustical signature have been 
observ'ed during conditions of light rainfall (Nystuen and Farmer, 1987). The surface 
wind effect will be described later. Second, the character of the underwater sound 
generated by rainfall appears to be influenced not only by the amount of rainfall, but 



by the sizes of the raindrops that strike the surface. This possibility was suggested by 



McGlolhin as one explanation for the inadequate empirical fit between rainfall rate 
and underwater spectral levels which he obtained for conditions of moderate to heavy 
rainfall. 

The purpose of this thesis is to predict the underwater sound spectrum observed 
at sea during conditions of moderate to heavy rainfall using laboratory measurements 
of the undervv'ater sound generated by large raindrops. Chapter II will explain the 
motivation for this work. In the first part of Chapter II, the mechanisms for under- 
water sound generation by raindrops will be discussed. Characteristics of raindrop 
size distributions will then be described, along with implications of drop size 
distribution variations on rainfall rate predictions using underw'ater sound measure- 
ments. The remainder of Chapter II will describe field ob.servations of underwater 
rain noise and relate these ob.servations to the raindrop sound generation mechanisms. 

In Chapter III. results of laborator}’ measurements of the impact noise of 
raindrops will be presented. The primary' purpose of these measurements was to 
determine how to treat the laborator>’ raindrop impact sound in terms of its farfield 
contribution to the underw ater rainfall acoustical noise. In Chapter IV, the acoustical 
energ) density spectrum associated with various sizes of large raindrops will be 
obtained from laboratory measurements of large raindrop impact and bubble sound. 
In Chapter V, the relationship between the laboratory measured spectra and the 
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rainfall spectra observed at sea will be established, and then used to predict the 
underwater sound spectrum due to rainfall. An inversion technique for determining 
raindrop size distributions from the measured rainfall sound spectrum is also 
described in Chapter V. The inversion technique is presented mainly to establish a 
formalism for future work. As will be seen, because of the lack of adequate informa- 
tion about the sea surface conditions, the rainfall drop size distributions, and the 
bubble densities during rainfall, an inversion was not attempted in this work. 
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II. BACKGROUND 



A. SOURCES OF UNDERWATER SOUND DUE TO RAINDROPS 

In order to determine the relationship between underw'ater acoustic levels and 
rainfall rate, the characteristics of the underwater sound produced by individual 
raindrops must first be understood. The tw o components of the underwater sound 
produced by drop splashes were first described by Franz (1959). The first source of 
raindrop sound is the impact itself; the second source is an oscillating bubble which 
sometimes forms. The underwater acoustical signal generated by a large (4.2 mm 
diameter) and particularh energetic drop is show n in Figure 2. The oscillating bubble 
is a highly efficient radiator of acoustical energy, generating a damped sinusoidal 
signal much longer in duration than the impact signal. For the case in Figure 2. the 
bubble formation is delayed by about 50 msec from the time of impact. 

Two di.stinct mechanisms have been identified for the formation of bubbles by 
natural!) occurring raindrops (Snyder, 1990). These are known as the Type I mech- 
anism (or regular entrainment) for small raindrops (0.8 to 1.1 mm diameter) and the 
Type II mechanism for large raindrops (greater than 2.2 mm diameter). The regions 
of bubble formation (in terms of drop size and impact velocity) for each of these two 
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Figure 2. Impact and bubble signal of a 4.2 mm diameter raindrop in artificial 
seawater. The time delay between impact and onset of the bubble is about 50 msec. 

mechanisms are shown in Figure 3. Note that for natural rainfall the drops strike the 

surface at terminal velocity, so only regions on the terminal velocity curve are 

applicable. 

The Type 1 mechanism occurs for drops 0.8 to 1.1 mm in diameter falling at 
terminal velocity, and has been studied at the Naval Postgraduate School (Kurgan, 
1989; Medwin et al., 1990) and elsewhere (Pumphrey et a]., 1989; Oguz and Prosper- 
etti, 1990). The Type 1 bubble results when the base of a conical splash crater is 
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Figure 3. Bubble formalion regions as a function of drop diameter (horizontal axis) 
and drop velocity (vertical axis) (Jacobus, 199 1). NCPA Ls the National Center for 
Physical Acoustics. 

pinched off (Longuet-Higgins, 1990); the resulting bubble has a resonance frequency 
of about 15 kHz. The Type I raindrops have been observed to produce bubbles 100% 
of the time when striking a smooth water surface at normal incidence, but the 
percentage of bubble formation decreases to 0% as the incidence angle increases to 25® 
(Kurgan, 1989).' 

The Type II mechanism occurs for drops larger than 2.2 mm in diameter. So far 
this mechanism has been studied only at the Naval Postgraduate School (Snyder, 
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1990; Jacobus, 1991) because of the availability of a 26 meter high drop tower (the 



larger Type II drops must be released from sufficient height to achieve terminal 
velocity). The sequence of events leading to Type II bubble formation are shown in 
Figure 4 (Snyder, 1989). First, a closed canopy of water forms above the splash 
cavity. Water continues to flow up the sides of the canopy. This convergence of water 
generates upw'ard and downward moving turbulent jets. The downward jet plunges 
through the bottom of the splash crater, entraining air that pinches off to form a 
resonating bubble. 

Unlike the Type I raindrops at normal incidence, not all drops larger than 2.2 
mm diameter produce bubbles when striking the surface at normal incidence and at 
terminal velocity. Also, the dependence of Type II bubble formation on incidence 
angle is unknown. However. Snyder (1990) observed that, in order for a bubble to 
form, the downward moving jet must be canted (as illustrated in Figure 4). Drops for 
w hich there was no cant in the downward moving jet were not observed to produce 
bubbles. This suggests that increasing the incidence angle may actually increase the 
probability of bubble formation, since oblique incidence may contribute to asymmetry 
in the resulting splash. 

The frequency spectrum of underwater sound generated by Type II raindrop 
bubbles was studied by Jacobus (1991). The spectral density curve of an individual 
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Figure 4. Drop sequence sketch (Snyder. 1990). 
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Type II raindrop bubble signal is shown in Figure 5. The figure shows the conlrib- 
ulion to the spectrum from the dominant bubble and from a secondary bubble which 
sometimes forms (more than one secondary bubble may be present). According to 
both Jacobus and Snyder, Type II dominant bubbles typically resonate at frequencies 
between 2 and 10 kHz. The frequency of a freely oscillating bubble in water is 
inversely proportional to the bubble radius (Qay and Medwin, 1977). The dominant 
Type II bubble frequency tends to decrease with increasing drop size (Figure 6), 
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Figure 5. Spectrum of the 1 m on-axis bubble signals for a 3.4 mm raindrop in 
artificial seawater combined with filtered seawater. The peaks in the spectrum are due 
to primary and secondary bubbles. 



10 



suggesting that larger raindrops lend to produce larger bubbles as expected. Jacobus 
also studied the variation of Type II raindrop bubble energy with drop temperature, 
surface temperature and salinity. He determined that Type II bubbles radiate more 
energy as the temperature difference between the drop and the surface increases. In 
fact, the average bubble energy for a 10° C temperature difference is almost twice that 
for a 0° C temperature difference. Also, the bubbles radiate much more energy in 
fresh water than in saline water (he reports a 45% lower bubble energy in water of 35 
ppl salinity as compared to fresh water). Both of these effects have significant 
implications for the prediction of rainfall rates at sea from acoustical measurements. 
The cau.se of the bubble energy variations with salinity and .surface-drop temperature 
difference has not yet been determined. 

B. CHARACTERISTICS OF RAINFALL DROP SIZE DISTRIBUTIONS 

Because the character of the underwater noise generated by individual raindrops 
is a function of the drop .size, it may also be true that the character of underwater 
rainfall sound is dependent on the number of drops of a given drop size that strike the 
surface. If underwater acoustical signals are to be u.sed for rainfall mea.surements, the 
raindrop size distribution must be known or inferred, as well as the relationship 
between the raindrop size distribution and rainfall rate. A summary of raindrop sizes 
and the corresponding mechanisms that generate underwater noise is presented in 
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Drop Diameter (mm) 



a Jacobus (1991) + Snyder (1 990) 

Figure 6. Type II dominant bubble frequency vs. drop diameter. 
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Table 1. The nomenclature shown in Table 1 for the various drop size ranges was 
suggested by Medwin et a!. (1992). 

TABLE 1. RAINDROP SIZES VS. SOURCES OF ACOUSTIC NOISE 



Nomenclature 


Equivalent Raindrop 
Diameter (mm) 


Source of Underu^ater 
Acoustic Noise 


Minuscule 


O 

O 

bo 


Impact Only 


Small (Type I) 


0.8 - 1.1 


Bubble and Impact 


Mid-Size 


1.1 - 2.2 


Impact Only 


Large (Type II) 


2.2 and larger 


Bubble and Impact 



The classic raindrop size distribution model is the Marshall and Palmer model 
( 194<S) given by; 

.V( D ) = ■ e 

where A' is the number of drops per unit volume (m^) per unit drop diameter 
increment (0. 1 nun). D is the equivalent spherical drop diameter (mm), and .V, (drops 
per nr' per 0.1 mm) and A (mm ') are statistical parameters that determine the 
character of the drop size distribution. The term "equivalent spherical diameter" is 
used here since, e.xcept for minuscule drops, raindrops are not in fact spherical; they 
are ellipsoidal, with the larger ones flattened at the base. Marshall and Palmer 

determined that the parameter A (the exponential slope) is given by A = 4.1 • R 
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where R is the rainfall rate in mm/hr. Marshall and Palmer also suggested that 
(the D = 0 intercept) is a constant equal to 800 drops per per 0.1 mm. Thus, the 
Marshall'Palmer distribution predicts that the number of very small drops present in 
rainfall is independent of the rainfall rate {N approaches as the drop size, D, goes 
to zero), and that fewer large drops are present for low rainfall rates than for high 
rainfall rates. The characteristics of the Marshall-Palmer distribution are illustrated in 
Figure 7, which shows the Marshall-Palmer distributions for rainfall rates of 15 and 
100 mm/hr. 

If the Marshall-Palmer drop size distribution were suitable for all rainfall 
conditions, then the rainfall rate could be determined by measuring the exponential 
slope of the drop size distribution alone, and then using the Marshall-Palmer 
Equation for zl to solve for R. However, actual raindrop size distributions can differ 
markedly from the Marshall-Palmer distribution. Cataneo and Stout (1968) observed 
drop size distributions on the Atlantic Coast which varied with synoptic conditions 
(cold frontal rains consisted of a larger number of small drops than warm frontal rains 
for similar rainfall rates). Waldvogel (1974) observed large jumps in the value of 
associated with mesoscale variations within a given precipitation field. For instance, 
A(, values for showers and thunderstorms (in areas of weak convection) tended to be 
large (on the order of 35,000) while values for steady, widespread rain (in areas 
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15 mm/hr 100 mm/hr 

Figure 7. Marshall-Palmer drop size distributions ,N {D ), for 15 mm/hr and 100 
mm/hr rainfall rates. 

with little or no convection) tended to be small (on the order of 4000). The value for 
A. the drop size distribution slope, has also been observed to depart from the 
expression given by Marshall and Palmer. Hodson (1986) observed that yl approaches 
a constant value between 2.1 and 2.3 mm ‘ for rainfall rates in excess of 25 mm/hr; 
above 25 mm/hr, an increase in rainfall rate is reflected by an increase in the value of 
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While the exponential distribution is accepted by most researchers to be an 
adequate approximation to actual raindrop size distributions, Ulbrich (1983) 
contends that drop size distributions can be better modeled as gamma distributions of 
the form; 

N{D) ^ D>‘ ■ (2) 

where the variable /i is a positive or negative value. Figure 8, obtained from Ulbrich 
(1983), compares the exponential distribution with the gamma distribution for both 
positive and negative values of /x. When /x is greater than zero, the distribution is 
concave downward, and the smallest and largest drops are fewer in number than for 
the exponential distribution (in fact, the number of drops goes to zero as D goes to 
zero for drops less than 1 mm diameter). When /x is less than zero, the distribution is 
concave upward, and the number of small and large drops is greater than for the 
comparable exponential distribution. Ulbrich also calculated the values of /x and 
for dilTerent environmental conditions by empirical analysis of rainfall data gathered 
from many different sources. He determined that the drop size distribution param- 
eters can be roughly categorized according to the type of precipitation. For instance, 
thunderstorms tend to have relatively large values of and positive /x, while is 
smaller and /x more variable for stratiform precipitation. 
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Figure 8. Gamma raindrop size distributions, N{D), for positive and negative values 
of ji. where D is the drop diameter (Ulbrich, 1983). 

WTiile many measurements of raindrop size distributions have been conducted in 
connection with the study of radar reflectivity from rainfall, very few concurrent 
measurements of rainfall drop size distributions and rainfall underwater acoustical 
signatures have been made. Because of this, it will later be necessary to assume a form 
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for the drop size distribution in order to calculate the predicted underwater acoustic 
signal from the laboratory measured acoustic signal of individual raindrops. The form 
of drop size distribution for any required computations will be the exponential 
Marshall-Palmer distribution. However, it must be kept in mind that the actual drop 
size distribution may differ significantly from the Marshall-Palmer distribution, which 
in turn may cause significant differences between the actual and predicted acoustical 
rainfall signature. Simultaneous field measurements of rainfall rate, rainfall drop size 
distribution, and rainfall acoustic signals are needed to remedy this limitation. 

The raindrop size distribution can be integrated to yield the rainfall rale 
according to the following equation; 



w here R is the rainfall rate in mm/hr. N{D) is the raindrop size distribution (drops 



in mm (Hodson. 1986). The terminal raindrop velocity is a function of raindrop 
diameter, and varies from approximately 4 m/s for raindrops of 1.0 mm equivalent 
spherical diameter to 9.1 m/s for raindrops of 5.0 mm equivalent spherical diameter 
(Pruppacher and Klelt, 1978). Snyder (1990) demonstrated that, for raindrops of size 
2.2 to 5 mm diameter, the terminal velocity can be approximated by the equation *» 
4.6 • D''^ (m/s). Thus, for raindrops in excess of 2.2 mm diameter, the rainfall rate 




( 3 ) 



per m’ per 0.1 mm) is the terminal drop velocity (m/s), and D is the drop diameter 



18 



integrand varies approximately as the 7/2 power of D. Since the contribution to the 
total rainfall rale increases rapidly with increasing drop size, the contribution of large 
raindrops to the total rainfall rate can be significant, even though the raindrop size 
distribution tends to decrease exponentially. Since large raindrops also generate a 
characteristic underw-ater acoustical signature (via the Type II bubble generating 
mechanism), they w'ill be emphasized in this work. 

C. FIELD OBSERVATIONS OF RAINFALL ACOUSTIC SIGNATURES 

1. Observations for Light Rainfall 

Much of the research on the ambient underwater noise generated by 
rainfall has focused on light rainfall conditions because of a consistently observed 
peak in the underwater sound spectrum at about 15 kHz when light rain is present. 
An example of the underwater sound spectrum for light rain is shown in Figure 9 
(Nystuen and Farmer. 1987). The 15 kHz frequency in the spectrum for light rainfall 
has been well correlated with the number of small (< 1.5 mm diameter) raindrops 
present (Nystuen. 1986). The source of the spectral peak at 15 kHz has been mainly 
attributed to the generation of Type 1 bubbles by raindrops of 0.8 to 1 .1 mm diameter, 
since the Type I bubbles have a characteristic resonance frequency of about 1 5 kHz 
(Pumphre) et al.. 1989; Medw'in eta].. 1990). 
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Figure 9. Undei^’aler sound speclnim during light rainfall showing characteristic 1 5 
kHz peak (Nystuen and Farmer, 1987) 

The underwater spectral levels in the vicinity of 15 kHz are not useful for 
predicting the total rainfall rate, since the magnitude and shape of the ambient rainfall 
noise peak at 15 kHz depends not only on the number of 0.8 to 1.1 mm drops present 
in the rainfall, but also on the wind speed and the surface roughness. This is 
illustrated in Figure 10 (Nystuen, 1992), which shows the ambient noise generated by 
0,6 mmdir rainfall for various w'ind speeds. As the w-ind speed increases from 0.6 m/s 
to 3.3 m/s, the spectral peak near 15 kHz becomes less pronounced. The reason for 
the observed change is that, as wind speed increases, a greater percentage of the small 
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frequency in klU 

Figure 10. Effect of wind on the 15 kHz spectral peak during light rain. As wind 
speed increases, the 15 kHz peak becomes weaker. Predicted spectra (based on the 
bubble formation percentages in Table 2) are superimposed (Nystuen, 1992). 

raindrops strike the surface at oblique incidence, and. as discussed earlier, the 

likelihood of Type I bubble formation for oblique incidence is very' small. In fact, the 

percentage of bubbles formed at the various wind speeds shown in Figure 10 was 

calculated by N\stuen (1992). The calculation assumed a form for the surface 

roughness based on the ambient wind speed; the results are presented in Table 2. 

Clearly, the presence of even light to moderate winds can have an overwhelming effect 

on the underwater sound generated by the Type I raindrop bubbles. 
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TABLE 2. THE PERCENTAGE OF DROPS CREATING BUBBLES 
FOR DIFFERENT WIND SPEEDS (NYSTUEN, 1992) 





Wind Speed (m/s) 


0 


0.6 


1.4 


2.6 


3.3 


Drop 

Size 

(mm) 


0.8 


100 


53 


12 


0.3 


0.02 


0.9 


100 


58 


15 


0.8 


0.1 


1.0 


100 


61 


18 


1.4 


0.3 



2. Observations for Moderate to Heavy Rainfall 

The characteristics of the underwater acoustic signature for moderate to 
heavy rainfall rates is much different from that for light rainfall. Based on measure- 
ments of rainfall underwater acoustical spectra at the Ocean Test Platform (OTP) in 
the Gulf of Mississippi, Tan (1990) di,scovered that the 15 kHz peak observed in the 
spectrum for light rainfall is no longer present during conditions of moderate to heavy 
rainfall. This is illustrated by the spectrum in Figure 1 1, obtained during conditions of 
heavy convective precipitation at the OTP (McGlothin, 1991). 

A positive correlation between the number of large raindrops (greater than 
2.2 mm diameter) and underwater spectral levels at lower frequencies (less than about 
10 kHz) has been recognized for sometime. Nystuen (1986) observed at Clinton Lake, 
Illinois that the underwater rain noise levels at frequencies less than 10 kHz (and at 
higher frequencies as well) increased as the number of large raindrops increased. 
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Figure II. Underwater sound spectrum during heavy precipitation (McGlothin, 
1991 ). . 



Scrimger ei al. ( 1987) also observ ed an increase the underw ater spectral levels at low 
frequencies w ith an increase in the number of large raindrops present. At the time of 
those studies, the mechanism for the increase in spectral levels at low frequencies w'ith 
an increase in the number of large drops was unknown. It is now known that the Type 
II bubble generation mechanism is the most likely source of the low frequency noise 
generated by large raindrops. 

Recent measurements of the correlation between rainfall rate and under- 
water spectral levels were performed by McGlothin (1991). The results are shown in 
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Figure 12. The correlation coefficients obtained by McGlothin are uniformly high 
over most of the frequency spectrum, except during events when the wind speeds are 
high (greater than 10 m/s). At the higher frequencies, the correlation coefficients 
between rainfall rate and sound level decrease for high wind speed. This is thought to 
be due to bubble clouds generated by breaking waves, a phenomenon which has been 
obsen'ed in the field (Farmer and Lemon, 1984), and studied in the laboratory 
(Medwin and Daniel, 1990). 

McGlothin's results for events with rainfall rates greater than 1 50 mm/hr 
are shown in Figure 13. For frequencies between about 2 kHz and 10 kHz, the 
correlation coefficient is uniformly high (greater than 0.8) for all the events shown, 
even during high wind conditions. Note that for higher frequencies, the correlation 
coefficient decreases for all events (although the decrease is most marked with high 
wind present). McGlothin suggests that this may be due to subsurface bubble clouds 
generated by the heavy rainfall itself, w'hich would also attenuate the high frequency 
rainfall noise as it propagated through the cloud. 

Given the strong correlation betw'een rainfall rates and underw'ater sound 
levels, especially below' 1 0 kHz, and a knowledge of the sound production mechanisms 
of individual raindrops, it should be possible to use measurements of the underwater 
sound generated by rain to infer the numbers of raindrops of various sizes striking the 
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surface. Assuming that all or mo.st of the rain noise that is < 10 kHz is generated by 
large (Type II) raindrops, a drop size distribution for the large drops can be estimated, 
and then extrapolated to smaller raindrop sizes. The drop size distribution obtained in 
this way can then be integrated to yield the total rainfall rate. This is the nature of the 
inverse problem of obtaining rainfall rate from the acou.stical signature of the rainfall. 
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Correlation CocfTicienl 




Figure 12. Correlation coefTicients between underwater spectral levels and rainfall 
rate for the rainfall events studied by McGlothin (1991). 
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Correlation Coefficient 




Figure 13. Correlation coefficients between underwater spectral level and rainfall 
rate v\here the rainfall rate exceeded 150 mm hr (McGlothin. 1991 ). 
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III. CHARACTERISTICS OF THE IMPACT SOUND 



A. BACKGROUND 

The impact sound of individual raindrops contributes to the total underwater 
noise generated by rainfall at sea. Figure 14 shows a typical impact acoustic signal 
obtained using a hydrophone 6 cm below the surface with no electronic filtering. The 
two components to the impact sound are a short duration impulse and a relatively 
long duration plateau. The corresponding frequency spectrum is shown in Figure 15. 
In the frequency domain, the impulse component shows up as broadband noise and 
the plateau shows up as a low frequency spike. Pumphrey (1991) points out that, due 
to the broadband nature of the impact spectrum, any measurements of impact sound 
will be affected by the impulse response of the measuring equipment. To demonstrate 
this, the impact signal shown in Figure 14 was filtered with a 1 to 30 kHz bandpass 
filter (as will be explained later, these are the equipment filter settings used when 
measuring the bubble signals). The results are shown in Figures 16 and 17. Figure 16 
shows the filtered impact signal, while Figure 17 shows the corresponding frequency 
spectrum. Note that the peak impulse pressure of the filtered signal is significantly less 
than that of the unfiltered signal, and that the plateau component of the filtered signal 
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Figure 14. Laboraiory measured underwater acoustical signal for a 4.6 mm diameter 
raindrop impact with no filtering. The hydrophone was at a depth of 5.5 cm. 



is diminished. A detailed empirical study of the impact sound of raindrops was 
performed by Pumphrey and Crum (1989). They determined that the impulse 
component of the impact tended to var>’ with range as 1/r, while the plateau compo- 
nent tended to fall off more rapidly v\ith range (approximately as 1/r). Pumphrey 
(1991) treats the plateau component as a near field acoustical effect, and uses this 
assumption to derive the corresponding far field pressure associated with the plateau. 
This approach contradicts the results of Nystuen (1986) who, based on numerical 
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Figure 15. Frequency spectrum of the unfiJlered impact signal showm in Figure 14. 

simulations of a raindrop impact, asserted that the plateau is a hydrod>mamic effect 
associated with fluid flow in the vicinity of the impact. 

Laboratory' measurements of the raindrop impact sound have been conducted 
with three objectives in mind. The first objective was to investigate some of the results 
obtained by other researchers pertaining to characteristics of the impact sound. The 
second objective was to resolve how' to best treat the plateau component of the impact 
noise in terms of its contribution to the far field acoustic pressure (if any). The third 
objective was to determine the average impact frequency spectra for the various sizes 



30 




Figure 16. Effeci of applying a 1 lo 30 kHz bandpass filter to the impact signal. The 
filter causes the signal to go negative for a short time after the initial impulse. 



of large raindrops as a basis for solution of the inverse problem (measuring rainfall 
rate from the unden^ater acoustical signature of large raindrops). The data collection 
and analysis techniques for obtaining the impact spectra were similar to those used to 
obtain the bubble spectra, and will be described later. 
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Figure 1 7. Frequency spectrum of the filtered impact signal. 

B. EXPERIMENTAL PROCEDURES 

1. Pressure vs. Range Dependence 

To determine the pressure versus range dependence of the impact and 
plateau components, the e.xperimental setup shown in Figure 18 was employed. A set 
of 25 raindrops of 4.6 mm diameter were launched from a height of 26 m in a 
ventilation shaft that serves as a drop tower. The 26 m drop height ensures that the 
large 4.6 mm drops strike the surface at terminal velocity (as does natural rainfall). At 
the base of the tower is a redwood tank, which is 1.5 m deep and 1.5 m in diameter. 
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The lank is lined with redwood wedges lo reduce reverberation from the tank walls. 
For this experiment, the lank was filled with fresh water. The temperature of the tank 
water was 30° C. The device used to generate the drops consisted of a calibrated 
pipette tip attached to an intravenous drop bottle, with a drop accuracy of ± 5 percent 
b\ volume. 

To measure the acoustic signal, two hydrophones were used. The first 
hydrophone was located at a depth of 5.5 cm. and consisted of a laboraton,' con- 
structed hydrophone with a lead zirconate acoustic element and a nominal sensitivity 
of -206 dB re V'^Pa. The second h\drophone was an LC-10 located directly below 
the first hydrophone and at a depth of 22 cm. The nominal sensitiviiv of the lower 
hydrophone w as - 1 9S dB re V |iPa. The h\ drophone signals were then amplified by a 
factor of 100 using HP465.^ amplifiers, and patched to Computerscope. a digital 
multichannel analyzer mounted in an IBM PC XT computer. Computerscope's 
maximum sampling rate is 1 MHz divided by the number of input channels used (up 
to 16 channels are available) The maximum rate of 500 kHz was used for sampling 
the signal from the two hy drophones. 

The difference in the arrival time of the impact signal at the two hydro- 
phones was used to calculate the range and angle from the point of impact to the 
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Figure 18. Schematic diagram of the experimental setup for the impact pressure vs. 
range measurement. 
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position of the hydrophone. The geometry of the problem is shown in Figure 19. The 



horizontal distance A' in cm is given by; 



A = 



h:-h:-c--ai- 






2 -c-At 



- h; 



(4) 



w here A i is the arrival time difTerence (sec), c the speed of sound in water (cm/sec). H, 
and H, the hydrophone depths (cm), and c-Ai the range difference from the impact 
to the tw o hydrophones. After solving for .V. the range.s R, and /?. and the angles 6, 



and 6 to the respective hy drophones were determined from the equations 



hr, 



and 6 = arctan(.\' H). 

2. Radiation Pattern Measurement 

The impact radiation pattern was investigated using the setup shown in 
Figure 20. For this e.xperiment. an array of 4 hydrophones yvas constructed. The 
hydrophones yvere attached to a semicircular ring of 8 cm radius, and positioned at 
angles of0°. 20'’. 40°. and 60° yvith the vertical. The hydrophone output signals w'ere 
again amplified b\ a factor of 100. and patched to Computerscope. yvhere the signals 
yvere sampled at a rate of 250 kHz. Raindrops of 4.6 mm diameter were dropped from 
a height of 3 m into an anechoically lined fresh yvater tank containing the hydrophone 
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Figure 19. Geometry of the impact pressure vs. range measurement. 

array. The reduced drop height was necessary to ensure that a sufficient number of 

drops landed directly over the center of the array. 



3. Peak Impact Pressure vs. Impact Velocity 

To investigate the relationship between impact velocity and impact 
acoustic pressure, another experiment was conducted using a setup similar to that for 
the radiation pattern measurement, but using only one hydrophone directly below the 
point of impact. The impact velocity of the drops was varied by adjusting the height 
of the dropper. The relationship between dropper height and impact velocity is given 
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Hgure 20. Schematic diagram of the experimental setup for impact radiation pattern 
measurement. 



37 



by the following equation; 



K, = v^- {\ - exp (^^)) (5) 

w'here v, is the impact velocity in m/s, Vj. is the terminal velocity in m/s (9. 1 m/s for 4.6 
mm diameter drops), g the gravitational constant in m/s\ and h the drop height in m 
(Pumphrey and Crum, 1989). 

C. RESULTS 

1 . Impact Radiation Pattern 

The impact radiation pattern was measured for two reasons. First, for the 
pressure vs. range experiment, the measured impact pressures needed to be corrected 
to the corresponding on-axis values. Second, the radiation patterns of the impulse and 
plateau components of the impact sound needed to be compared to determine if there 
were any differences. For a simple acoustic source in water located near the surface, 
the radiation pattern may be expected to be that of a dipole, due to the presence of a 
virtual image of the .source above the surface of the w-ater and 180° out of phase with 
the source. This effect is illustrated in Figure 21 . The far field pressure radiated by an 
acoustic dipole is given by: 

P{r,ej) - - • cos d • (6) 

r 
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Figure 21. Negative image of a simple acoustic source near a pressure release 
boundarv (the air-water interface). Resulting pressure field is dipolar in nature. 

where D is the dipole source strength (Pa • m). r is the range, 6 is the angle between 

the dipole a.xis and the range vector to the far field point, and the expionential term 

gives the phase variation with distance and time. Thus, the acoustic pressure radiated 

by a compact dipole source region of this nature should vary as cos 0, where 6 is the 

angle with the vertical below the source. If the radiation pattern does not exhibit a 

cos 6 dependence, the source region may be of a different nature (such as monopole 

or quadropole). or the pressure field may be of nonacoustic origin. 

Results of the radiation pattern measurements are shown in Figures 22 and 

2."^ for the impulse and plateau components respectively. The figures show the 

normalized pressure vs. angle for three different impacts, along with the theoretical 
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dipole radiation pattern. The pressures have been normalized to the on axis (zero 
degree) value. The radiation patterns of both the impulse and plateau tend to conform 
to the cos 6 dependence of a dipole. This can also be seen in Figure 24, which shows 
the pressure time series obtained at the 0“ and 40“ hydrophones for a single impact. 
The two time series have been superimposed, and the 40“ hydrophone signal has been 
corrected by dividing by the cosine of 40 degrees. The magnitudes of the super- 
imposed signals for both the impulse and the plateau are nearly identical. 




Figure 22. Polar radiation pattern of the normalized impulse pressure for three 
impacts. The dipole axis is horizontal. 
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Figure 23. Polar radiation pattern of the normalized plateau pressure for three 
impacts. The dipole axis is horizontal. 

V^Tiile the impulse and plateau components of the impact sound both 
appear to have a dipole radiation pattern, it is not necessarily true that the plateau 
pressure is of acoustical origin. If the plateau pressure is due to incompressible 
(hydrodynamic) fluid flow associated v^iih the impact, then the cos 6 dependence may 
result from the fluid flow below the impact being predominantly in the vertical 
direction, since the projection of the fluid velocities into directions away from the 
vertical will also var>' as cos 6. A numerical study of raindrop impacts conducted by 
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Nystuen (1986) supports the latter contention. Results discussed later refer to the 
pressure versus velocity finding. 

2. Impact Pressure vs. Range Variation 

The impact pressure time series obtained at the upper and lo>ver hydro- 
phones for the pressure vs. range experiment are shown in Figure 25. The upper 
hydrophone was 5 cm below the surface, while the lower hydrophone was 22 cm below 
the surface. The impact shown landed directly above the hydrophones, so any 




Figure 24. Impact signal in freshwater at the 40® hydrophone (dashed), corrected for 
cos(40®), superimposed on the impact signal at the 0° hydrophone (solid). Range to 
the hydrophones was 8 cm. 
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difference between the two h\drophone signals can only be attributed to the differ- 
ences in range and in hydrophone sensitivities. Note that the relative magnitude of the 
impulse and plateau components is very different for the two hydrophones. At the 
lower hydrophone, the plateau component is barely discernible, while at the upper 
hydrophone, the plateau component is about 13 as large as the impulse component. 
This suggests that the impulse and plateau components do in fact have different range 
dependencies. 
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Figure 25. Pressure signals at the upper (top) and lower (bottom) hydrophones for 
an impact directly overhead. 
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Plots of axial pressure vs. range for both the plateau and peak impulse 
pressure are shown in Figure 26 for 4.6 mm diameter raindrops at terminal velocity. 
The pressure values were calculated using the following equation: 



byd 



axial 



M • G-cos d 



(7) 



where is the output voltage, M is the hydrophone sensitivity (V/Pa), and G is the 



amplifier voltage gain (100). The cos 6 term in the denominator is used to correct for 
the dipole angular variation of the pressure fields. The plateau pressures were 
recorded 200 iisec after the peak impulse pressure. The line fitted through the impulse 
data has a slope corresponding to a Mr range dependence, while the line fitted 
through the plateau data has a slope corresponding to a 1/r^ range dependence. 

While there is much scatter in the data, the difference in range dependence 
between the two components of the impact sound can readily be seen from the 
increasing spread between the impact and plateau pressures with range. The impulse 
pressure varies approximately as 1/r, while the plateau pressure varies approximately 
as l/r‘. The data points for the deeper hydrophone, however, tend to fall below the 
theoretical curves. One reason for this may be differences in response of the two 
hydrophones to the very rapid pressure change associated with the impulse (on the 
order of 4 to 8 iisec). 
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Figure 26. On-axis pressure vs. range for the impulse and plateau components (4.6 
mm drops). The hand fitted curves show a Mr dependence (impulse) and 1/r^ 
dependence (plateau). 
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The rapid decrease of plateau pressure with range indicates that the plateau 
is a near field effect of either acoustic or hydrodynamic origin. The Mr variation of 
the much larger impulse pressure, though, agrees with the expected far field pressure 
variation of an acoustic dipole. Therefore, the contribution of the plateau to the 
measured signal at a distant hydrophone will be negligible when compared with the far 
field contribution of the impulse. 

3. Dependence of Impact Pressure on Raindrop Velocity 

Franz (1954) asserted that the impact pressure of a raindrop would be 
similar to that of a rigid sphere .striking a water surface, for which the initial impact 
pressure at a range r below the sphere would be given by: 



P = ( 



p cos 6 






rc 



( 8 ) 



w here r is the impact velocity, r is the range, a is the radius of the sphere, c is the 
sound speed, p is the fluid density, and is the depth of penetration of the sphere. 
The peak impact pressure w ould occur at time zero, and w ould be proportional to 
A numerical analysis of raindrop impacts (Nystuen, 1986), how'ever, revealed that the 
rigid sphere model is oversimplified. One flaw is the sphericity assumption; actual 
large raindrops more closely resemble oblate spheroids that are flattened at the base. 
Nystuen's numerical results indicate that the drop shape can have a significant effect 
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on ihe underwater pressure generated by a drop, as shown in Figure 27 (the peak 
impulse pressure predicted numerically is about three times as large for the flattened 
drop as it is for the spherical drop). Nystuen also observed that the initial impact 
pressure in the immediate vicinity of the impact varied as v rather than as for a 
given drop shape. 

Laboratory measurements of initial impact acoustic pressure (impulse 
pressure) as a function of terminal velocity tend to agree with Franz's theory. The 
drop shape, however, cannot be varied independently of the velocit> for the laborator>' 
measurements. Pumphrey (1989) empirically determined the peak impulse pressure of 
the impact to var) as the 2.7 power of velocit> . His experiments were performed u.sing 
drops of 2.52 and .7.8 mm diameter. Figure 28 shows the results of the peak impulse 
pressure vs. impact velocity experiment using 4.6 mm diameter drops. TTte slope of the 
line fitted to the data corresponds to pressure vary ing as the .7.5 pow'er of the velocity. 
The reason the value obtained here is larger than Pumphrey’s may be due to increased 
flattening of the large 4.6 mm diameter drop with increasing velocity, since the 
flattening is more pronounced for larger drops (Pruppacher and Pitter. 1971). The 
reason why Nystuen's numerical model of raindrop impacts fails to correctly predict 
the impact pres.sure's dependence on impact velocity is uncertain, but it may be due to 
the fact that Nystuen's numerical simulations of impact pressure w'ere for field points 
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in the immediate vicinity of the impact (at distances on the order of mm) rather than in 
the far field of the impact source. The characteristics of the proximal sound field for 
the impacts are expected to differ significantly from the far field characteristics. 

The dependence of the peak impulse pressure on the size and velocity of a 
raindrop has important implications for the relative contribution of different size 
raindrops to the total underwater sound generated by rainfall. Assuming a d 
dependence of peak pressure on drop diameter and terminal velocity, a 5.0 mm 
diameter drop will have an impulse peak pressure 54 times (or 35 dB greater than) that 
of a 1 .0 mm diameter drop. Thus, despite their smaller numbers, the impact sound 
radiated by the larger raindrops (in excess of 2 mm diameter) can make a significant 
contribution to the total rainfall impact noise. 
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B) flAttimad drop 



Figure 27. Numerical prediction of the effect of drop shape on the on-axis pressure 
directly beneath the point of impact for a 3.0 mm diameter raindrop . A realistic 
flattened drop shape is shown in the lower figure (Nystuen, 1986). 
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Figure 28. On-axLs peak impulse pressure vs. impact velocity for 4.6 mm diameter 
raindrops in freshwater. The hydrophone depth was 1.5 cm. 
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IV. ENERGY SPECTRA FOR LARGE RAINDROPS IN SALTWATER 



A. PURPOSE 

As mentioned earlier, based on the work of Snyder (1990) and Jacobus (1991) 
there is reason to believe that much of the underwater noise generated by rainfall 
during moderate to heavy rainfall conditions is due to bubble and impact sound 
generated b\ large raindrops. Therefore, additional measurements have been made of 
the energy spectral levels of the underwater sound due to individual raindrops, over 
the size range of drops which produce Type II bubbles (2.2 mm diameter and larger). 
Once the underwater sound energy generated by individual raindrops is known as a 
f unction of drop size, the contribution to the underwater sound from individual drops 
can be summed to obtain a predicted acoustic signature for a given rainfall rate. 

Measurement of the underwater sound energy spectra for various sizes of large 
raindrops was first accomplished by Jacobus (1991). One of the average spectra 
obtained by Jacobus for 4.2 mm diameter drops is shown in Figure 29. which shows 
the impact contribution to the intensity spectrum as well as the total spectrum (bubble 
and impact). Note that the impact spectrum is virtually insignificant when compared 
to the spectrum that includes the bubble noise. 
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Figure 29. Impact and bubble components of the intensity spectrum for 4.2 mm 
diameter raindrops, as measured by Jacobus (1991). The impact spectrum is very close 
to the X axis. 



The spectra obtained by Jacobus, however, are insufTicient for predicting 
underwater acoustic levels for rainfall. A careful review of Jacobus’ work revealed 
that he applied a Hamming window to his measurements of the bubble pressure 
signals prior to computing the intensity spectra of the bubbles. While his purpose for 
windowing his pressure signals was reasonable (to minimize leakage of the bubble 
energy to sidelobes in the frequency domain), he failed to account for the effect that 
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using a Hamming window would have in reducing the total measured intensity of the 
signal. For stationary random processes, the coherent proces.sing lo.ss (Harris, 1976) 
due to use of a Hamming window can be compensated for, but the bubble pressure 
.signal is transient in nature, and the usual coherent loss factor cannot be directly 
applied. 

To a.ssess the effect of applying a Hamming window to a Type II bubble pressure 
signal, an idealized signal was used. The idealized bubble pressure signal is given by: 

/>( t) = • cos(27tF/) 

where is the bubble resonance frequency, and a is the amplitude decay rate. For the 
test case. F was taken to be 2000 Hz and a u as taken to be (5 msec)'*, a value close to 
the theoretical decay rate for a 2000 Hz bubble in water. The analytic .solution for the 
Fourier Transform of the pressure signal is given by: 

nD ( 10 ) 

{2-kjT + a)‘+ {2ti To)' 

and so it is possible to compute the theoretical intensity spectrum of the ideal infinite 
duration bubble signal for the sake of comparison. A discrete finite form of the ideal 
pre.ssure .signal was generated using MATLAB (a mathematics program), with a 
simulated sampling rate of 250 kHz and a total record length of 1 6.4 ms (4096 points). 
The .simulated pres.sure signal is shown in Figure 30. Figure 31 shows the same signal 
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Figure 30. Theoretical bubble signal with a Rectangular window. 

after a Hamming window has been applied. The energy spectra of both the unwind- 
owed (or rather, the Rectangular windowed) signal and of the Hamming windowed 
signal were calculated u.sing the Discrete Fourier Transform, and compared to the 
theoretical intensity spectrum based on the analytic expression for the Fourier 
Transform. The result is shown in Figure 32, which shows the spectral peak at 2000 
Hz due to the bubble. The spectrum obtained using the rectangular window is very 
similar, both in shape and magnitude, to the theoretical spectrum. The spectrum 
obtained using the Hamming window, however, differs markedly both in shape and 
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Figure 31. Pressure signal after applying a Hamming window. 

magnitude from the theoretical bubble spectrum. In fact, the ratio of the total sound 
intensit) using the Hamming window to the total sound intensity using the Rectangu- 
lar window is about 0.13. Clearly, the effect of using a Hamming window makes 
.lacobus's results for the bubble spectra quantitatively unreliable. 

.Mso. the number of bubble signatures collected by Jacobus for each raindrop 
size was relatively small from a .statistical .standpoint (he collected about 30 samples 
for each drop size). To obtain a more reliable estimate of the average spectrum for 
each raindrop size, the goal here was to obtain 100 bubble signatures for each of six 
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Frequency (Hz) 



Figure 32. Comparison of the theoretical Fourier spectrum of the pressure signal 
(solid line) with the Rectangular windowed DFT spectrum (crosses) and the Hamming 
v\ indowed DFT spectrum (dashed). 

large drop sizes. The individual bubble spectra were then averaged and smoothed to 
obtain representative spectra for various drop size ranges. A summary of the drop 
diameters, size ranges, and number of bubble samples obtained is presented in Table 3. 
The drops of each diameter listed in Table 3 are taken to be characteristic of the 
corresponding range of drop sizes. Over 90 .samples were obtained for all but the 
.smalle.st drop size (the low bubble production rate, and deflection by drafts in the 
ventilation shaft, made collection of the 2.5 mm drop bubble signals difficult). 
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Since the impact spectra are much smaller in magnitude than the bubble spectra 
for each drop size, only 30 impact samples were collected to obtain average impact 
spectra for three drop diameters (2.5 mm. 3.4 mm. and 4.2 mm). The impact spectra 
for the remaining drop sizes were obtained by interpolating the measured spectra (the 
method used will be described later). 



TABLE 3. NUMBER OF BUBBLE SAMPLES FOR EACH DROP SIZE 



Drop Diameter (mm) 


Drop Size Range 
(mm) 


Number of Samples 


2.5 


2.2 - 2.8 


46 


3.0 


2.8 - 3.2 


99 


3.4 


3.2 - 3.6 


98 


3.8 


3.6 - 4.0 


100 


4.2 


4.0 - 4.4 


100 


4.6 


4.4 - 4.8 


99 



B. EXPERIMENTAL PROCEDURES 

1 . Setup for Bubble Sound Measurements 

The equipment setup used to measure the Type II bubble signal is shown in 
Figure 33. All the drops w'ere released from a height of about 26 m using the drop 
tower and redwood tank arrangement described in the previous chapter. For these 
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experiments the tank was filled with saline water. Initially, artificial sea salt was added 
to the tank to simulate ocean water. Later, filtered seawater was obtained from the 
Monterey Bay Aquarium and added to the water already in the tank. 



IV Dropper (Larger Drops) 

or Pipette (Smaller Drops) 

V 



26 m 



I 




O i 

O I 



PC/AT ^ 

with Computerscope (A/D) 



tthaco Preamp 




Krohn-Hrte Filter 



1.5 m i 

o • 

Hydrophone 

AnechoicTank ' 




Figure 33. Equipment setup for bubble signal measurements. 

For the drops in the 3.4 mm to 4.6 mm size range, calibrated pipette tips 
attached to an I-\^ bottle were used to generate the drops, with an accuracy of ± 5% by 
volume. For the other two drop sizes (2.5 mm and 3.0 mm), the drops were too small 
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lo fall on their own. At first, the smaller drop sizes w'ere generated by attaching a 
pipette tip to a loudspeaker driven by a function generator operating in "burst" mode. 
When the function generator triggered, the drop would be forced off the pip^ette by the 
vibration of the loudspeaker. The accuracy of this method was determined lo be ± 
20%, which was unsuitable when compared to the accuracy obtained for the larger 
drop sizes. The bubble signal measurements for the smallest tw'o drop sizes were 
therefore repeated using an Eppendorf digital pipette, with an accuracy of ± 1% by 
volume. The disadvantage of u.sing the pipette is that the drops had to be generated 
manually, requiring tw'o people for data collection. 

A single hydrophone was used to measure the underwater bubble noise. 
The hydrophone contains two coaxial 1/8 inch cylindrical barium lilanate elements . 
and has a nominal sensitivity of -91 .5 dB re V/Pa. The response of the hydrophone is 
flat(± 3 dB) over a frequency range of 5 to 300 kHz (Snyder, 1990). The hydrophone 
w'as located at a depth z of 9 cm for the 4.2 mm and 4.6 mm drops, and at a depth z 
of 6.0 to 6.5 cm for the smaller drops. The distance h on the surface from the point 
above the hydrophone to the point of drop impact w'as measured with a ruler to an 
estimated accuracy of ± 0.5 cm. A metal grid was placed over the surface lo help 
determine where the drops landed; the grid squares were 5 cm by 5 cm. 
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An Ithaco preamplifier was used to amplify and filler the bubble signals. 
A gain of 2000 was used for the larger drop sizes and 5000 for the smaller drop sizes. 
The Ithaco was set up as a band pass filter with a pass band of 1 to 30 kHz. The filter 
settings were necessary because of electrical interference in the building at frequencies 
outside the pass band range. The output of the Ithaco was then passed through a 
single Krohn-Hite band pass filter with the same filter settings to further attenuate 
interfering noise. The resonance frequencies of the bubbles were in the pass band 
range, and the filler settings did not interfere with the bubble signal measurements. 

The data was collected using the Computerscope digital analyzer, now 
mounted in a PC/AT computer. The sampling rate for the bubble signals was 250 
kHz. The length of the bubble records extracted for processing and analysis was 16.0 
ms (or 4000 data points for each bubble). A 16 ms record length was necessary to 
recover most of the energy (over 95%) of each bubble signal measured. 

Measurements were also made of the drop temperature, surface tempera- 
ture, salinity, and surface tension. The temperature measurements were made with 
mercury thermometers, w'ith an accuracy of ± 0.5° C. The salinity was measured 
using a salinometer accurate to 0.05 ppl. Surface tension was measured with a 
capillary' tube, which was accurate to ± 5 dyne/cm. 
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As mentioned in the introduction. Jacobus (1991) discovered that the 



energy radiated by the Type II bubbles seemed to depend on the salinity and on the 
surface-drop temperature difference. A summary of the present temperature 
measurements is given in Table 4. While the tank temperature remained fairly 
constant, the fluid temperature of the drops was difficult to control, since the data 
collection process could take several hours. Subsequently, the drop - surface 
temperature difference varied between 2 C° and 5.5 C°. This could result in a ± 0.8 
dB variation in spectral energy, based on Jacobus's data for fresh water. 



TABLE 4. SUMMARY OF DROP AND SURFACE TEMPERATURES 



Drop Diameter 
(mm) 


Drop Tempera- 
ture (°C) 


Surface Temper- 
ature (°C) 


Temperature 

Difference 

(C°) 


2.5 


28 


23 


5 


?.o 


21 


23 


2 


3.4 


24 


21 


3 


3.8 


24 


21 


3 


4.2 


26 


23 


3 


4.6 


17 


22.5 


5.5 
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A summar>' of the salinity and surface tension measurements is shown in 



Table 5. The salinity values ranged from about 32.3 ppt to 37.3 ppt. The work by 
Jacobus suggests that thus could cause a ± 3% (or ± 0. 1 dB) variation in the energy 
levels. The tank saluuty increased after the filtered seawater was added, apparently 
due to salt which had collected on the walls of the tank where the level was lower. 
Evaporation of some of the tank water may also have contributed to increasing 
salinity levels. In general, the surface tension in the tank increased with time due to 
repeated filtering of the tank water and the addition of pure filtered seawater after 
collection of the 4.2 mm raindrop data. Also, the water in the tank is discolored by 
the leaching of tannic acid from the redwood. Overall, the surface tension of the tank 
water was 10 to 20 dynes/cm lower than the surface tension of pure water at 18° C (73 
dynes/cm). The effect this had on the bubble sound or the generation of bubbles is 
unknown. 

2. Setup for Impact Sound Measurements 

While the setup used to obtain the impact spectra was similar to the setup 
for the bubble measurements, there were some important differences. First, the 
impact signals were not band pass filtered. Filtering was minimized to prevent altering 
the impact signal. The Ithaco preamplifier was used to amplify the impact signals 
using a gain of 500, and had a high pass filter setting of 0.03 Hz to remove dc bias 
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TABLE 5. SUMMARY OF SURFACE TENSION AND SALINITY 



Drop Diameter (mm) 


Salinity (ppt) 


Surface Tension (dyne/cm) 


2.5 


37.3 


62 


3.0 


37.2 


62 


3.4 


34.9 


57 


3.8 


34.9 


57 


4.2 


34.1 


57 


4.6 


32.3 


52 



from the input signal to Computerscope. The Krohn-Hite band pass filter was not 
used at all. The hydrophone was placed at a depth of 10 cm for the impact measure- 
ments. 

The impact signals were collected using Computerscope at the maximum 
sampling rate of 1 MHz. The largest possible sampling rate w'as used to prevent 
undersampling the impulse component of the impact sound. The noise from 
interference was much more pronounced than for the bubble measurements, and the 
impact signals were much shorter in duration. Therefore, shorter length records were 
extracted for the impact signals than for the bubble signals. The extracted record 
lengths were 500 psec for the 2.5 mm impacts, 1 msec for the 3.4 mm impacts, and 2 
msec for the 4.2 mm impacts. 
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C. DATA ANALYSIS 



1 . Analysis of the Bubble Signals 

a. Conversion of Hydrophone Voltage to On-axis Pressure 

The hydrophone voltage levels recorded by Computerscope were 
corrected to the on axis pressures at 1 m. To accomplish this, the recorded voltages 
were first corrected for the amplifier gain and the hydrophone sensitivity using: 



Pbyd 



V 

G • M 



( 11 ) 



where is the pressure at the hydrophone (Pa), V is the recorded voltage, G is the 
amplifier gain, and M is the hydrophone sensitivity (V/Pa). The hydrophone response 
was assumed to be omnidirectional. 

Second, a correction must be applied to account for divergence {Mr 
pressure variation), and a dipole radiation pattern (cos 6 pressure variation). The 
pertinent geometr>’ is shown in Figure 34. The dipole radiation pattern for bubble 
oscillations has been confirmed by Kurgan (1989) for Type I (small raindrop) bubbles. 
The necessar>- correction factor is given by; 



on O.WS = - 

100 ‘Z 



( 12 ) 
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where is the uncorrecled pressure, z is the hydrophone depth (in cm), b is the 
distance from the point of impact to the point above the hydrophone (in cm), and 
Piwonixis ^he on-axis pressure at one meter. 




Figure 34. Geometry for range and angle correction. 

For some of the bubble signal measurements, the hydrophone w-as in 

the near field of the dipole. To correct for this, a third correction factor had to be 
applied, and is given by (Medwin and Beaky, 1989): 



P unhid ” P tm on xxis 



N 



1 



(13) 



where is the corrected far field pressure, r is the range (m), and k is the wave 

number (m *). which can be computed from the resonance frequency of the bubble {k 
- 2-Kflc, where c is the sound speed in m/s and f the bubble frequency in Hz). 

b. Calculation of the Acoustic Energy Spectra 

A 4096 point Fast Fourier Transform of each bubble signal time series 
was computed using MATLAB to obtain the bubble frequency spectra. The FFT 
algorithm employed by MATLAB is given by: 

\_1 -i7nnk 

A'(L+1) = ^ x{n-^\)e ^ 

/ 7-0 

where .v {n) represents the time series values. X {k) represents the frequency series 
values, and N is the number of points in the series. For this case, the x {n) are 
pressure values in Pa; the units of the spectral values X{k) are Pa per bin. 

To compute the energy density spectrum in Pa^* s/Hz from the FFT 
values, the following equation is used: 



E{f) - - 



At ' LML)P 
N ■ df 



(15) 



w here E {f) is the energy spectral level (in units of Pa^* s/Hz) at the frequency f, X 
the Discrete Fourier Transform value with index k, and N is the number of points. 
The frequency f is given by • df as k varies from 0 to N !2. The frequency 
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resolution (or bin width) is given by df - MNAU where 4/ is the sampling interval. 
For the 250 kHz sampling rale used to collect the bubble data, the sampling interval 
was 4 ^sec, giving a frequency resolution of 61 Hz for a 4096 point FFT. To resolve 
any confusion between Equation 15 and other equations for the power or energy 
spectrum of a signal, a derivation of Equation 15 is provided in Appendix A. 

c. A veraging of the Bubble Spectra 

The individual bubble sound spectra will contain peaks at the 
resonance frequencies of the bubbles. More than one peak may be present if the event 
results in secondary bubbles as well as a dominant bubble. Both the magnitude and 
the frequency of the peak(s) can var>' from one spectrum to the next. It is therefore 
neces.sar> to statistically average the bubble spectra in some fashion to obtain a 
representative spectrum for a given drop size. 

To accomplish this, two different approaches have been employed 
here. The first approach, following the method of Jacobus (1991). is to ensemble 
average the individual bubble spectra for each raindrop size. After the spectra are 
ensemble averaged, a 1 kHz wide moving filter is applied to smooth the ensemble 
spectrum. 

An alternative approach for obtaining a smooth spectrum was to fit a 
smooth curve to the probability distribution of dominant bubble frequencies for each 
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drop size. The bubble frequency distribution curv'e n {f), where n (/’) represents the 
fraction of bubbles with resonance frequencies within a given frequency interval 
centered on the frequency 4 is then multiplied by the total spectral energy (in units of 
Pa^- s) to obtain a smooth energy density spectrum for each drop size. The total 
spectral energy for each drop size was obtained by integrating the average spectra 
obtained previously. 

The smoothing method just described is performed using only the 
dominant bubble frequencies, since the dominant bubbles seem to account for most of 
the raindrop bubble energy (Jacobus, 1991). To show this, the average spectra were 
computed using the dominant bubble contribution only and compared to the spectra 
containing both dominant and secondar) bubble energy. To obtain only the dom- 
inant bubble contribution to the average spectra, a simulated bandpass filter centered 
on the dominant bubble frequency was applied to the individual bubble spectra prior 
to averaging. The filter frequency width was six times the theoretical 3 dB bandwidth 
for an oscillating bubble at a given frequency. This bandwidth was selected because it 
allowed for recovery’ of approximately 90% of the dominant bubble energy. 
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d. Bubble Total Energy Calculation 



To investigate the variation of dominant bubble energy with bubble 



frequency, the dominant bubble energies were calculated using the following equation 
(Kurgan, 1989): 



where E is the total bubble energy (in Joules), r is the range, pc is the acoustic 
impedance, and the summation is over the square of the axial pressure values for the 
bubble signal time series. This equation accounts for the dipole nature of the bubble 
sound, and holds for k ■ r> 3, where k is the wavenumber. A derivation of Equation 
16 is given in Appendix B. The squared axial pressure summation in the time domain 
is related to the axial energy density spectrum defined in Equation 15 by; 



Equation 16 can therefore be expressed in terms of the axial energy density spectrum 
as; 




(16) 



5: £-(/)J/'= Ep.,(0^ A/ 



(17) 




3 pc- /.(» 



(18) 
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where the range r is 1 m. The integration is performed over a frequency interval equal 
to six times the 3 dB bandwidth of the bubble signals. This bandwidth was sufficient 
for recovering 90% of the dominant bubble energy. 

2. Analysis of the Impact Signals 

Similar analysis techniques were used to obtain the impact spectra, with 
some important differences. For converting the measured voltages to pressures, no 
near field correction factor was applied, since the hydrophone was placed farther from 
the surface than for the bubble measurements, and the impact noise is broadband. 
Also, the record lengths used for the spectral analysis were shorter (from 500 jisec to 
2 msec). The contribution to the impact spectra from the plateau component of the 
impact noise, which results in a peak in the spectrum at low frequencies (less than 
about 1000 Hz), was rejected. As stated earlier, the plateau noise does not appear to 
contribute to the far field noise signature of the impact. 

D. RESULTS 

1 . Bubble Formation Percentages 

A summary of the percentage of large raindrops observed to produce 
bubbles in salt water is given in Table 6. The percentages are in general agreement 
with the percentages obtained by Snyder (1990) for fresh water. In general, the 
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percentage oflarge raindrops that produce bubbles tends to decrease with decreasing 



drop size. 

TABLE 6. BUBBLE PRODUCTION PERCENTAGE 



Drop Diameter (mm) 


Bubble Percentage 


Fresh 


Salt 


2.5 


20 


25 


3.0 


37 


37 


3.4 


50 


50 


3.8 


63 


50 


4.2 


62 


50 


4.6 


55 


62 



or the fraction of raindrops \\ hich produce bubbles, a still smaller fraction 
produce both dominant and secondar) bubbles. The percentage of secondar)’ bubble 
formation for the drops obsen ed to produce bubbles is given in Table 7. Overall, the 
percentage of .secondary bubble formation is between 1 0 and 20%. The relatively high 
percentage ob.served for the 4.6 mm drops could be because they are the most 
energetic. As noted by Jacobus (1991 ), the .secondary bubbles were generally smaller 
in magnitude and higher in frequency than the dominant bubbles. 
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TABLE 7. PERCENTAGE OF BUBBLE PRODUCING DROPS THAT 
INCLUDE SECONDARY BUBBLES IN SALT WATER 



Drop Diameter (mm) 


Percentage 


2.5 


8 


3.0 


9 


3.4 


19 


3.8 


16 


4.2 


12 


4.6 


35 



2. Average Bubble Spectra 

The spectra obtained from ensemble averaging and smoothing the 
individual bubble spectra are shown in Figures 35 through 40. Even after smoothing, 
the spectra have multiple peaks. Jacobus ascribed the secondar>' peaks in the spectra 
to the contribution of secondar>’ bubbles and to a phenomenon he referred to as 
"wobbles". The wobbles were gated sinusoidal signals that appeared before, during, or 
after the bubble signal. No wobbles were observed in the bubble signals collected for 
this work; in fact, the wobbles may have been a tank scattering artifact which was 
present only prior to installation of the redwood wedges on the walls of the tank. 
Moreover, the likelihood that the secondary bubbles would produce peaks in the 
spectra of the same order of magnitude as the dominant peak is inconsistent with the 



72 



spectral Level (Pa ''2 * s/Hz) 
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relatively small percentage of secondary bubble formation, and the lower energy of the 
secondary bubbles as compared to the dominant bubbles. 




Figure 35. Average bubble spectrum, 4.6 mm diameter raindrops. 
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Figure 36. Average bubble spectrum, 4.2 mm diameter raindrops. 
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Figure 37. Average bubble spectrum. 3.8 mm diameter raindrops. 
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Figure 38. Average bubble speclrum, 3.4 mm diameter raindrops. 
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Figure 39. Average bubble spectrum, 3.0 mm diameter raindrops. 



77 



spectral Level (Pa'^2 * s/Hz) 
(Times 10E-9) 




0 ^ ^ ^ ' 

0 2 4 6 8 10 12 14 16 18 20 



Frequency, kHz 



Figure 40. Average bubble spectrum, 2.5 mm diameter raindrops. 
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To investigate the notion that the spectral peaks are due mainly to the 
dominant bubble signals, average spectra were obtained which included only the 
energy of the dominant bubble. The resulting spectra tend to have the same overall 
shape as before, although the spectral levels are somewhat smaller (about 10 to 20%) 
over most of the frequency range. An example is shown in Figure 41 for the 4.6 mm 
diameter drops. Similar results were obtained for the other drop sizes. 




Figure 41. Comparison of the 4.6 mm raindrop bubble spectrum (solid) with the 
spectrum of the dominant bubble energy alone (dashed). 



79 



The broad shape of the spectra can be mostly attributed to the dominant 
bubbles. However, highly energetic individual bubbles may be responsible for local 
peaks in the spectra, since the spectral peak for an individual bubble can be relatively 
narrow (the 3 dB bandwidth for a 2000 Hz bubble is on the order of 100 Hz). An 
ensemble average of as many as 1 00 bubble spectra may be insufficient for obtaining a 
smooth spectrum. All of the average bubble spectra have been plotted in Figure 42 to 
illustrate their relative shapes and magnitude. 
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Figure 42. Bubble spectra for the six large raindrop diameters studied. 
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3. Bubble Frequency Distributions 

The dominant bubble frequencies for each drop size were computed. The 
mean, the median, and the standard deviation of the dominant bubble frequencies are 
given in Table 8. As can be seen, the dominant bubbles tend to increase in frequency 
with decreasing drop size. Also, the variation in bubble frequencies becomes quite 
large for the smallest drop sizes (2.5 mm and 3.0 mm). Both of these observations 
agree with the trends observed in the average bubble spectra. 



TABLE 8. DOMINANT BUBBLE STATISTICS 



Drop Diameter 
(mm) 


Mean Frequency 
(Hz) 


Median 

Frequency (Hz) 


Standard 
Deviation (Hz) 


2.5 


10,023 


9125 


5597 


3.0 


9874 


8484 


6442 


3.4 


4905 


4242 


2305 


3.8 


4514 


3998 


2139 


4.2 


3716 


3235 


2240 


4.6 


3327 


2808 


1999 



Histograms showing the percentage of dominant bubbles vs. frequency are 
shown in Figures 43 through 48. The frequency bins in these figures are 244 Hz wide. 
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The smooth, filled curves shown in the figures are gamma probability dislribulions, 
given by: 

g{f) ? — /'>0 

P‘r(a) 

where /is the frequency, a and B are parameters that define the shape of the distribu- 
tions. and r(a) is the gamma function evaluated at a (eg., see Walpole and Meyers, 
1989). The gamma distribution was selected to model the dominant bubble frequency 
distributions because the distributions appear to be skewed right (towards higher 
frequencies). The skewness could not have been modeled adequately using a normal 
probability distribution unless higher moments (skewness and kurtosis) had been 
introduced. Most of the filled curves are in reasonable agreement with the actual 
dislribulions. The fits are poorest for the smallest drop sizes (2.5 and 3.0 mm 
diameter). For the 2.5 mm drops, this may be in part due to the smaller sample size 
(46 vice 98 to 100 for the other drop sizes). Since the average energy for the smaller 
drops is less than for the larger drops, the fit may not be as critical for these sizes. 

The characteristics of the bubble frequency distributions as a function of 
drop size ma> explain some of the observ'ed spectral characteristics. As drop size 
decreases, the dominant bubble frequency tends to increase (eg., the smaller drops 
tend to produce smaller, higher frequency bubbles). The standard deviation of 
dominant bubble frequencies also lends to increase with decreasing drop size. In fact. 
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Figure 43. Dominant bubble frequency distribution, 4.6 mm diameter drops. The 
dashed line is the gamma distribution of Equation 19 with a = 4.2>\, B= 626. 
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Figure 44. Dominant bubble frequency distribution, 4.2 mm diameter drops. The 
dashed line is the gamma distribution of Equation 19 with a = 10.9. 303. 
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Figure 45. Dominant bubble frequency distribution. 3.8 mm diameter drops. The 
dashed line is the gamma distribution of Equation 19 with a - 1 1.3, /?= 372. 
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Figure 46. Dominant bubble frequency distribution, 3.4 mm diameter drops. The 
dashed line is the gamma distribution of Equation 19 with a - 10.2, B = 502. 
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Figure 47. Dominant bubble frequency distribution, 3.0 mm diameter drops. The 
dashed line is the gamma distribution of Equation 19 with a = 2.25, B = 4000. 
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Figure 48. Dominant bubble frequency distribution, 2.5 mm diameter drops. The 
dashed line is the gamma distribution of Equation 19 with a - 5.98, D- 1840. 
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the probability distributions for the smallest drop sizes (2.5 and 3.0 mm diameter) are 
rather uniform, which agrees with the observation that the spectra for the smallest 
drops are rather broad. In the next section, the fitted distribution curves will be used 
to obtain smoothed average spectra for the various drop size ranges. The smoothed 
spectra will then serve as the basis for calculating the predicted underwater acoustic 
spectra due to rainfall at sea for a given rainfall rate. 

4. Filled Bubble Speclra 

The fitted distributions for the dominant bubble frequencies were u.sed as 
the basis for smoothing the bubble energy spectra. The average bubble spectrum for 
each drop size was integrated, and the resulting value was multiplied by the fitted 
bubble frequency distribution to yield a smooth spectrum. The results are shown in 
Figures 49 through 54. \\here the smoothed spectrum is compared with the average 
spectrum for each drop size. 

The smoothing technique used here does have some flaws. For instance, 
theoretically the energy radiated by individual bubbles depends on the bubble 
frequency. The damping constant for higher frequency bubbles is higher than for 
lower frequency bubbles, varying from 0.03 for 2 kHz bubbles to 0.06 for 20 kHz 
bubbles (Clay and Medwin, 1977). Thus, for the same initial amplitude, a high 
frequency bubble signal will decay more rapidly and radiate less energy than for a low 
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Figure 49. Average spectrum (solid) vs. smoothed spectrum (dashed), 4.6 mm drops. 
The integrated spectral energy is 2.13 • 10'* Pa'-s. 
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Figure 50. Average .spectrum (solid) vs. smoothed spectrum (dashed), 4.2 mm drops. 
The integrated spectral energy is 3.22 • lO"* Pa"-s. 
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Figure 51 . Average spectrum (solid) vs. smoothed spectrum (dashed), 3.8 mm drops. 
The integrated spectral energy is 1.1 1 •10'^ Pa‘*s. 
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Figure 52. Average spectrum (solid) vs. smoothed spectrum (dashed), 3.4 mm drops. 
The integrated spectral energy is 1.22 • lO"* Pa^* s. 
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Figure 53. A\ erage spectrum (solid) vs. smoothed spectrum (dashed), 3.0 mm drops. 
The integrated spectral energy is 1.14 • 10“* Pa‘-s. 
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Figure 54. Average spectrum (solid) vs. smoothed spectrum (dashed), 2.5 mm drops. 
The integrated spectral energy is 2.29 • 10"^ Pa'*s. 
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frequency bubble signal. For instance, a theoretical 2 kHz bubble with the same initial 



amplitude as a 10 kHz bubble will radiate 8 times as much energy. To investigate this 
possibility, the dominant bubble energies for each drop size were calculated per 
Equation 16 and then averaged together in 2 kHz wide frequency intervals. Plots of 
the average dominant bubble energy vs. frequency are shown in the upper half of 
Figures 55 through 60 for each drop size. The number of dominant bubbles in each 
frequency inter\'al is shown in the lower half of these figures. No energy vs. frequency 
trend can be readily discerned for most of the drop sizes, perhaps because of the 
relatively small number of drops in some frequency intervals. It may be possible that 
the higher frequency bubbles tend to have larger initial amplitudes than the lower 
frequency bubbles. For the 3.4 mm drops, however, the bubble energies do appear to 
decrease with increasing frequency. 

Some of the average bubble energies in the range of 1 2 to 1 8 kHz appear to 
be unusually high for the 4.6 mm, 4.2 mm, and 2.5 mm diameter drops. Some of the 
bubbles in this frequency range may actually have been secondary raindrop bubbles 
generated by droplets expelled during the raindrop splash. 

Another factor to consider is the variation of the bubble signal bandwidths 
with frequency. The theoretical bandwidth for a bubble oscillating at a frequency 
is given by: 
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( 20 ) 



where is the 3 dB bandwidth and is the damping constant (Clay and Medwin, 
1977). A 2000 Hz bubble, for instance, has a theoretical bandwidth of about 60 Hz in 
the frequency domain, while the bandwidth for a 20 kHz bubble is about 1200 Hz. 
Thus, not all the bubble energy for the higher frequency bubbles will be confined to a 
single frequency bin (for the fitted distributions, the bins were 244 Hz wide). 

Despite these Haws, the smoothed spectra are in reasonable agreement with 
the average spectra, except for the 2.5 mm drops. However, as can be seen from the 
energy vs. frequency plot for this drop size (Figure 60), it appears that the bubbles in 
the 16 to 18 kHz range were much more energetic (by a factor of more than 10) than 
the bubbles in the adjacent frequency bins. The peak in the average spectrum at about 
16 kHz for this drop size is probably due to these highly energetic bubbles, and may 
not have been prominent if a larger number of samples had been used. As mentioned 
earlier, some of the high frequency bubbles may actually have been secondary bubbles 
formed by splash droplets. 

5. Average Impact Spectra and Total Fitted Spectra 

The average impact spectra for the 2.5, 3.4, and 4.2 mm drops are shown in 
Figure 61. The broadband character of the impact sound is evident. Figure 62 is a 
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Figure 55. Average dominant bubble energy vs. frequency for the 4.6 mm drops 
(upper curve). Lower curve is a histogram of the number of bubbles per frequency 
bin. TTie energy value at 9 kHz may be anomalous. 
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Figure 56. Average dominant bubble energy vs. frequency for the 4.2 mm drops 
(upper curve). Lower curve is a histogram of the number of bubbles per frequency 
bin. 
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Figure 57. Average dominant bubble energy vs. frequency for the 3.8 mm drops 
(upper cur\^e). Lov^er curv'e is a histogram of the number of bubbles per frequency 
bin. 
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Figure 58. Average dominant bubble energy vs. frequency for the 3.4 mm drops 
(upper curve). Lower curve is a histogram of the number of bubbles per frequency 
bin. 
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Figure 59. Average dominant bubble energy vs. frequency for the 3.0 mm drops 
(upper curve). Lower curve is a histogram of the number of bubbles per frequency 
bin. 
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Figure 60. Average dominant bubble energy vs. frequency for the 2.5 mm drops 
(upper curve). Lower curve is a histogram of the number of bubbles per frequency 
bin. 
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Figure 61. Impact energy spectra for the 2.5 mm, 3.4 mm, and 4.2 mm diameter 
drops. 
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plot of the impact signal energy (in pJ) versus drop diameter. The energy of the 
impact signals varies as the 4th power of the drop diameter (d) over this size range. 
This agrees with the result obtained by Kurgan (1989) and Jacobus (1991) that the 
radiated impact acoustic energy is proportional to the kinetic energy of the drops, 
given by Va mv~ (the mass m of a raindrop is proportional while the terminal 
velocity for large raindrops is proportional to d '^’ ). The curve shown in Figure 62 was 
interpolated to obtain energy values for the other 3 drop sizes studied. The ratios of 
the impact energies were then multiplied by the measured spectra to obtain interpo- 
lated spectra for the other drop .sizes. 

All the impact spectra were then smoothed (using a 1 kHz moving average) 
and added to the smoothed bubble spectra. A summarv of the total impact and fitted 
bubble spectra for each drop size is .shown in Figure 63. The impact noise has some 
effect on the magnitude of the resulting spectra at the higher frequencies. 

A summar)' of the f>eak spectral levels vs. drop size is presented in Table 9. 
The peak spectral levels for the smallest drops are at least an order of magnitude less 
than for the largest drops. The smaller drop sizes have less kinetic energy, and so less 
energy is available to drive the bubbles into oscillation. In general, the following 
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Impact Energy (pJ) 




Drop Diameter (mm) 



Figure 62. Impact signal acoustic energ\ (pJ) vs. drop diameter (mm). The curve is a 
least squares fit. 



changes take place in the spectra with decreasing drop size: 



• The spectral levels decrease. 

• The spectra shift to higher frequencies. 

• The spectra become broader. 

The spectra shown in Figure 63 will be used in calculating the predicted underwater 
spectral levels at sea due to rainfall. 
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TABLE 9. PEAK SPECTRAL LEVELS VS. DROP SIZE 



Drop Diameter 
(mm) 


Peak Spectral 
Level (Pa^-s/Hz) 


2.5 


6.09 • lO '" 


3.0 


4.32 ■ 10'" 


3.4 


1,59 • lO ' 


3.8 


1.88 • 10'' 


4.2 


6.88 • lO' 


4.6 


4.80 • 10 ' 
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Figure 63. Total fitted energy density spectrum for each large raindrop size. 
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V. INVERTING THE FITTED RAINDROP SPECTRA 



A. INVERSE THEORY 

1 . Terms and Definitions 

Inverse techniques are of tremendous importance in a number of scientific 
fields, including ocean acoustical tomography and seismology. In simplest terms, the 
general linear inverse problem consists of finding the desired solution to a set of /n 
equations in n unknowns. The inverse method described here follows the approach 
used by Wiggins (1972), who outlined an inverse technique for studying surface waves 
in the Earth. 

Let RS{f) represent the rainfall under\\'ater acoustic spectrum measured at 
sea, and DR Did) the number of drops of a given drop size d falling per square meter 
per second. Also, let the function A (/Tc/) represent the contribution to the underwater 
acoustic spectrum at the frequency /from raindrops of a given diameter d, such that: 

RSi f) - Y. -DR Did) (21 ) 

d 

where the summation is over all drop sizes. An expression for A if,d ) will be derived 
later in determining the predicted undenv'ater spectrum due to rainfall from the 
individual raindrop spectra. 
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Because the raindrop acoustic signatures have been measured only for 
distinct drop sizes, and because the measured spectral levels are also discrete, the 
function A (fM ) must be expressed in discrete form. A { f,d ) then becomes a matrix, 
where the matrix entry A^- gives the contribution to the underwater sound spectrum 
from a drop size range Ad centered on drop diameter d- , to the frequency bin Af 
centered on frequency In a similar fashion, the functions RS{f) and DRD{d) can 
be represented as vectors, where RS^ is the spectral value of a frequency bin centered 
at /^, and DRDj represents the drop rate density (number of drops per m^ per sec) for 
drops of a given size range centered on diameter d^ . In what follows, matrices will be 
referred to using a capitalized bold letter, while vectors will be referred to using an 
uncapitalized bold letter. The frequency index / varies from 1 to m and the drop size 
index / \ aries from 1 to n. 

Using the matrix notation, the discrete form of Equation 21 is: 

D 

Y^A,;drd^- rs,, i-\,m (22) 

1 

which can be characterized as a set of m linear equations in the n unknown variables 
DRD^ . The coefficients of the variables are given by A-. , and the solution of each 



no 



equation is rs^ . The set of equations can be expressed in terms of matrix algebra as: 



A drd- IS 

where A is an m by ;? matrix, rs is the w by 1 solution vector containing the 
measured unden\’ater spectral values at discrete frequencies and drd is the vector 
containing the unknown drop rate densities for the discrete drop sizes d^ . 

The A matrix represents the contributions of each drop size range to the 
underwater spectral values in each frequency bin, and is therefore related to the 
physical mechanisms of underwater sound generation by raindrops. Later, a form of 
the A matrix will be derived using the individual raindrop spectra obtained in the 
previous chapter. The A matrix values, however, may also be determined statistically. 
A .statistical method that could be used to obtain the A matrix coefficients will be 
discussed later when making suggestions for future work. 

2. The Singular Value Decomposition Method 

Singular Value Decomposition (SVD) is the method used by Wiggins 
(1972) to solve the linear system A 'drd = rs for the unknown vector drd However, 
unlike Wiggins problem, the linear system for our problem is overdetermined; that is, 
the number of unknowns (given by /? ) is less than the number of equations (given by 
m ). Put differently, the number of frequency bins in the acoustic spectrum exceeds the 



111 



number of drop sizes studied. Unless at least m - /? of the equations are linearly 
dependent, an overdetermined linear system will have no exact solution for the 
unknown vector drd In this case, the linear system is said to be inconsistent (Anton 
and Rorres, 1987). In contrast, for an underdetermined system (where the number of 
unknowns exceeds the number of equations) an infinite number of solutions exist. 

Even though the linear system for our problem is overdetermined, the 
Singular Value Decomposition method can still be used to obtain a solution. While an 
exact solution to the system of equations may not be possible, the SVD method will 
still yield a best fit solution to the system of equations in a least squares sense. That is, 
if the solution drd (obtained by Singular Value Decomposition) is multiplied by A to 
yield the vector rs\ then the summation S(r.s^ - rs'y will be minimized. 

The first step in the SVD method is to decompose the A matrix into the 
product of three different matrices, such that: 

A ^ U‘A- ^-4) 

The matrix U is an m by k matrix that contains the eigenvectors of length m 
associated with the columns of the matrix A. The matrix is a diagonal matrix, 
where the diagonal elements are the eigenvalues of the matrix A. The number of 
eigenvalues is Ar, where k ^ n, and n is the number of unknowns in the set of linear 
equations. The matrix is the transpose of the n hy k matrix that contains the 
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eigenvectors of length n associated with the rows of A Techniques for obtaining the 
eigenvectors and eigenvalues of a matrix A are presented in Anton and Rorres (1987). 
Singular Value Decomposition can be accomplished with the MATLAB program by 
using the function SVD{A ). 

Once the matrix A has been decomposed, it can be substituted into 
Equation 23. Letting ‘drd = drd\ Equation 23 becomes; 



UAdrd' = rs 

Following Wiggins, the vector drd* will be referred to as the parameter correction 
vector. Since the columns of the matrix U consist of eigenvectors that are orthonor- 
mal. the product • U yields the Identity matrix / The product of A'' -A (where 
A is the eigenvalue matrix and A ' is its inverse) is also the identity matrix / by 
definition. The inverse of the matrix A is obtained by simply inverting the diagonal 
elements of A Multiplying both sides of Equation 25 by A '^ * £/^and applying the 
rules of matrix algebra leads to: 



drd' > A '-U^ rs 



(26) 
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An estimate for the unknown vector drd can then be obtained from the 



solution for drd' from: 

drd “ V’dixT 

The above expression is only approximately true, because the product V * does 
not necessarily yield an identity matrix (the eigenvectors are in the columns of V and 
in the rows of and matrix multiplication involves multiplying the rows of the first 
matrix by the columns of the second matrix). The product K* is referred to as the 
resolution matrix R. The more the R matrix resembles an actual identity matrix (with 
diagonal elements equal to 1 and off diagonal elements equal to zero), the better the 
solution drd is resolved. 

3. EfTecl of Noise on the Inverse Technique 

A solution for drd' can also be obtained using fewer than the maximum 
number of eigenvalues. To truncate the number of eigenvectors and eigenvalues to a 
value k. all the columns of U and V with indices greater than k are eliminated, as 
well as all columns and rows of A with indices greater than k. The effect of truncating 
the decomposition matrices is to reduce the goodness of fit of the calculated under- 
water spectrum rs’= A - drd as compared to the measured spectrum rs. 
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If the linear system described by Equation 23 were a perfect system, that is, 
if there were no errors in either the vector /s or in the matrix A, then the best solution 
for drd would be the solution obtained using the maximum number of eigenvalues, 
since this solution would yield the closest fit between the predicted and measured 
spectral levels. However, errors in the linear system (regardless of their source) will 
affect the results of inversion. The cumulative errors in the linear system relating the 
underwater acoustic levels to the drop rate density of drops striking the surface will be 
referred to as the system noise, and wall be represented as a vector e. The noise can be 
due to either random or biased errors. Pos.sible sources of noise include; 

• Measurement errors. 

• Modelling errors. 

• Effects not accounted for (wind noise, for example). 



The equation for the linear system (Equation 23) can be rewritten so as to 
include the cumulative noise effects; 



A' drd + e = rs 

Sohing for A ’drd Equation 28 becomes; 

A- drd - rs - e 



(28) 



(29) 
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Thus, the solution drd obtained using the SVD method in fact yields a best fit to rs-e 
rather than the best fit to rs alone. Significant difTerences between rs and rs-e will 
lead to errors in the inverse solution drd It may be true that a less perfect fit to the 
spectrum, obtained using fewer than the maximum number of eigenvalues and 
eigenvectors, will actually yield a better solution for the vector drd This possibility 
must be considered when applying the inversion method. 

In the following section, a form for the A matrix will be presented. Then, 
using a drop rate density vector computed from the Marshall-Palmer raindrop distri- 
bution, the underwater spectral levels for a given rainfall rate will be calculated, and 
compared to spectra measured at sea for similar rainfall rate conditions. If the 
predicted and the measured underwater acoustic spectra are in reasonable agreement, 
then an inversion technique based on the derived matrix A should be possible. 

B. SOLUTION OF THE FORWARD PROBLEM 
1 . Calculations 

Jacobus ( 1991 ) derived an expression for the underwater acoustic spectrum 
due to rainfall measured at an upward looking hydrophone at some depth below the 
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surface. The expression is duplicated here in modified form; 



- Jin X DRD{d) X E{f,d) x sinV(/) x 1 m' (30) 

d 

where DRD{d) is the drop rate density (in units of drops/m^ s) for each drop size, 
E{f.d) is the energy density spectrum (in units of Pa^- s/Hz) for each drop size, and 
ij/ {f) is the beamwidth of the receiving hydrophone. TTie summation is over all drop 
sizes d. A derivation of Equation 30 is presented in Appendix C. The rainfall 
spectrum can then be expressed in units of dB re pPa^/Hz using the equation; 

RSLiE) - 10 • log[/?5(/)] +120 

Comparing Equation 30 with Equation 21 from the previous section, the 
function A {f,d) that relates drop sizes to spectral levels is determined to be; 



A{Rd) - 7t • sin*^(/) • E{Ed) • 1 m^ 
The matrix (discrete) form of A is given by; 



(32) 



= Ti • 1 m^ • sinV(/) • E 

Each column of the matrix E contains the average energy spectrum (in Pa^- s/Hz) for 
a given raindrop size. 

Energy spectra have been measured for six large raindrop sizes, where the 
fitted spectra contain 82 frequency bins between 1.2 kHz and 21 kHz. A trial E 
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matrix was constructed with dimensions i = 1 to 82 and j - \ io 6. The E matrix 
was then substituted into Equation 33 to obtain a trial A matrix. The A matrix was 
then multiplied by an estimate of the drop rate density vector for a given rainfall 

rate to yield a prediction for the underwater acoustic spectrum due to rainfall. 

The computed spectrum was then compared to actual rainfall noise spectra 
measured at sea for the same rainfall rate. The at-sea spectra were determined by Tan 
(1990) and McGlothin (1991) from data obtained at the Ocean Test Platform (OTP) 
facility in the Gulf of Mexico. The location of the OTP is shown in Figure 64. 

The OTP hydrophone is an ITC-3001 (circular piston) transducer with a 
beamwidth of 28 degrees at 17.5 kHz, and is mounted in a concrete foundation in 15 m 
of water. The directivity function for a circular piston transducer is given by 

= 2 -IJ, ( V )l/ V , where v = ka - sin 6 and J, is the 0th Order Bessel Function (eg.. 

see Ziomek, 1985). In the expression for v, ka is the product of the w'aven umber and 
the piston radius, and 6 is the angle betw'een the field point and the acoustic axis of 
the piston. The 3 dB beamwidth is given by a value of v - 2.2, and so the beamw'idth 
of the OTP hydrophone can be estimated as a function of frequency. The result is 
show'n in Figure 65. The estimated OTP beamwidth w'as used in computing the A 
matrix values. For frequencies below' 8 kHz, the beam pattern of the OTP hydro- 
phone is omnidirectional. 
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Figure 64. Location of the Ocean Test Platform. 
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Figure 65. OTP hydrophone beamwdlh vs. frequency. 

Estimates of the drop rate density vector were obtained for two rainfall 
rates, 1 5 mm/hr and 1 00 mm/hr. The estimates were computed using the Marshall- 
Palmer drop size distribution described in Chapter II. The values N(d) (number of 
drops per m^ per 0.1 mm diameter increment) from the Marshall-Palmer distribution 
were converted to drop rate density values using the equation; 

DRDid) - (^) 
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where d is the drop diameter (mm), v, Ls the drop terminal velocity (m/s), and Ad 
the drop size range (in tenths of mm) centered on diameter d. The results are shown in 
Table 10 for the two rainfall rates. 

2. Results 

The rainfall spectra obtained by multiplying the A matrix by the estimated 
drop rate density vectors are shown in Figure 66. It should be noted that the effects of 
small and midsize raindrops were not included in obtaining the predicted levels. The 
peak spectral level for the 100 mm/hr rainfall is predicted to be 74.5 dB (re ^iPa^/Hz). 
The peak level for the 15 mm^hr rainfall is predicted to be 62 dB. Both of the 
predicted spectra have a negative slope of 10 to 11 dB per octave for frequencies 
greater than 8 kHz. All estimates assume that measurements are made in the free field 
(eg., away from surfaces). 

The predicted underwater acoustic spectra are only in approximate agree- 
ment with the rainfall spectra measured at sea. Figure 67 shows the predicted 15 
mm/hr rainfall spectrum, along with rainfall sound spectra measured at the Ocean 
Test Platform for the same rainfall rate. The ocean measurements have been 
referenced to free field values by correcting for reflection from the rigid concrete pad 
surrounding the OTP hydrophone. This is done by subtracting 6 dB from the 
measured spectral levels for frequencies greater than 1 kHz. The basis for this 
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Table 10. DROP RATE DENSITIES VS. RAINDROP SIZE 





Drop Rate Density (drops/m^s) 


Mean 
Drop Size 


Size 

Range 


1 5 mm/hr Rainfall 
Rate 


100 mm/hr Rainfall 
Rate 


2.5 


2.2-2.8 


109 


729 


3.0 


2.8-3.2 


24.9 


241 


3.4 


3.2-3.6 


10.2 


135 


3.8 


3.6-4.0 


4.13 


74.2 


4.2 


4.0-4.4 


1.67 


40.6 


4.6 


4.4-4.8 


0.67 


22.7 



correction will be described later. The wind speeds during the various rainfall events 
are shown in the upper right comer of the figure. Wind speeds were moderate for all 
three cases. The spectral levels for two of the rainfall events are above the predicted 
spectral levels by 5 to 10 dB for frequencies greater than 8 kHz and by 0 to 5 dB for 
frequencies between 5 kHz and 8 kHz. As discussed later, the measured levels below 
about 5 kHz are judged to be unreliable due to probable errors in the OTP hydro- 
phone calibration. 

The spectral levels for the third rainfall event shown in Figure 67 are much 
lower than for the other two rainfall events, and are less than the predicted spectral 
levels by as much as 1 5 dB for frequencies between 4 and 1 2 kHz. It is speculated that 
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Figiire 66. Predicted rainfall spectra for rainfall rates of 1 5 mm/hr and 100 mm/hr. 

the third spectrum was measured during a period when the number of large raindrops 
present was significantly less than for the other two cases. 

A smaller difference between the predicted and the measured spectral levels 
at the higher frequencies was observed for the 100 mmltr rainfall rate. Measured 
rainfall spectra for three events are compared with the predicted 1 00 mm/hr rainfall 
spectrum in Figure 68. For the two events with moderate wind speeds (6 to 7 m/s), the 
measured spectra are 3 to 9 dB greater than the predicted spectrum for frequencies 
greater than 1 0 kHz, and are within 6 dB of the predicted spectrum for frequencies 



123 



between 5 kHz and 10 kHz. Once again, the measured spectral levels below about 5 
kHz are judged to be unreliable. 

The third measured spectrum in Figure 68 was for a rainfall event where 
the wind speed was high. The spectral levels for frequencies greater than about 10 kHz 
are significantly less for this third event than for the other two events. As mentioned 
in Chapter II, the attenuation of the rainfall sound for frequencies greater than about 
10 kHz in the presence of high winds is probably due to subsurface bubble plumes 
generated by breaking waves. 

As mentioned earlier, the small and midsize drops were not included in the 
computation of the predicted rainfall spectrum. The characteristic 15 kHz peak 
associated with the small (Type I) drops has not been observ'ed at the OTP during 
moderate to heavy rainfall conditions, indicating diminished production of Type I 
bubbles. Also the small drop bubble noise is dependent on wind speed. For these 
reasons, small drops were excluded from the rainfall spectrum computations. 

To estimate the contribution to the spectrum from the midsize drops 
(which only generate impact noise), the relationship between drop diameter and 
impact energy was extrapolated to four of the midsize drop diameters. The drop 
diameters used were 1.15 mm, 1.45 ram, 1.75 mm, and 2.05 mm. The extrapolated 
midsize impact energies were then divided by the impact energy for the 2.5 mm drops. 
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(ZH/Z ap) 

Figure 67. Predicted vs. measured rainfall spectra, 15 mm/hr rainfaU rate. WS is the 
wind speed. 
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(ZH/Z ap) 13^91 

Figure 68. Predicted vs. measured rainfall spectra, 1 00 mm/hr rainfall rate. WS is the 
wind speed. 
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Figure 69. EfTect of midsize raindrops on the predicted 15 mm/hr rainfall spectrum. 
The dashed curve includes the contribution of midsize drops. 

and the ratio multiplied by the 2.5 mm raindrop impact spectrum. The midsize impact 
spectra were used to augment the A matrix, which was then multiplied by the drop 
rate density vector which included the midsize drops. The result is shown in Figure 69 
for the 1 5 mm/hr rainfall rate and Figure 70 for the 1 00 mm/hr rainfall rate. Including 
the estimated eff^t of the midsize impacts raised the spectral levels by as much as 2 dB 
for the 15 mm/hr rainfall rate and 1 dB for the 100 mm/hr rainfall rate. Apparently, 
the effect of the midsize raindrop impacts on the underwater spectral levels cannot be 
ignored. 
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Figure 70. EfTect of midsize raindrops on the predicted 1 00 mm/hr rainfall spectrum. 
The dashed curve includes the contribution of midsize drops. 

3. E>iscussion 

Recall that the measured spectral levels for most of the rainfall events were 
larger than the predicted spectral levels for frequencies greater than about 1 0 kHz. The 
reason(s) for these differences must be determined. The factors that may be responsi- 
ble for this can be divided into two categories; those that cause the measured spectral 
levels to be too high and those that cause the predicted spectral levels to be too low. 
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One of the factors which could result in measured spectral levels that are 



too large is the scattering effect of the OTP hydrophone mounting. The OTP hydro- 
phone is mounted flush in a concrete block. TTie radius of the hydrophone’s receiving 
surface is about 10 cm, while the block’s maximum dimension on its upper surface is 
about 2 feet or 60 cm. The presence of the block has a tremendous effect on the 
pressure measured at the hydrophone. When the wavelength of incident sound is 
significantly less than the dimensions of the block, the block will behave as a rigid 
baffle with respect to the hydrophone. For a simple receiver in the vicinity of an 
infinite baffle, the measured intensities are increased by a factor of 6 dB. The 
frequency corresponding to a half-wavelength of 0.6 m (the maximum dimension of 
the block) is about 1 kHz. However, the beam pattern of the hydrophone becomes 
directional above about 8 kHz. The overall effect that the hydrophone mounting has 
on the pressure measured at the hydrophone is therefore a complicated function of the 
frequency. The correction of 6 dB to the measured spectral levels for frequencies 
greater than 1 kHz is an approximation to this. An attempt must be made to more 
accurately quantify the effect of the hydrophone mounting on the measured levels. 

Another factor which could result in higher than actual measured levels are 
hydrophone calibration errors. These errors will be most pronounced below 5 kHz, 
the frequency below' which no hydrophone sensitivity values were provided by the 
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hydrophone manufacturer. An in-silu calibration of the ITC-3001 hydrophone has 
been conducted, and the calibration data is currently being analyzed with the intent of 
obtaining reliable hydrophone sensitivity values for the frequency range of interest. 

A number of factors may be responsible for making the predicted sound 
levels too low. The first is the effect of ocean surface roughness. The effect of surface 
roughness on the sound generated by individual large raindrops has not yet been 
studied. It may be possible, for instance, that the bubble formation percentages 
increase in the presence of surface roughness, and may even result in Type II bubble 
formation extending to drop sizes less than 2.2 mm diameter. Work is currently in 
progress here at NPS to measure the influence that surface roughness has on large 
raindrop underwater sound. 

Another factor that would make the predicted levels too low is the presence 
of a greater number of large raindrops in the raindrop size distributions than given by 
by the Marshall-Palmer distribution. Many of the rainfall events studied by Tan 
(1990) and McGlothin (1991) occurred during periods of convective rainfall. As 
discussed in Chapter II, the characteristics of raindrop size distributions can differ 
markedly from the Marshall-Palmer distribution and in fact depend on the type of 
precipitation. A drop size measuring device needs to be acquired in order to measure 
drop size distributions and underwater spectral levels concurrently. 
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A third factor that may be responsible for making the predicted spectral 
levels too low is the failure to account for the acoustic propagation conditions at the 
Ocean Test Platform. The OTP hydrophone Ls in a shallow water environment in 
which reverberation (eg., multipath arrival) must be considered. The expression for 
the A matrix was derived assuming free field conditions. It Ls recommended that open 
ocean measurements of the rainfall underwater sound spectrum be obtained, since the 
spectra measured in the open ocean should be in closer agreement with the spectral 
levels that are predicted using the free field derivation. For now. an attempt must be 
made to take into account the propagation conditions at the OTP. 

Another factor which may make the predicted levels too low are the effects 
of water property variations on the energy radiated by individual raindrops. Since the 
laborator>' measurements were performed using seawater, the effect of differences 
between ocean salinity and laboratory' tank salinity should be small. Jacobus (1991) 
empirically determined that the large raindrop acoustic energy is proportional to the 
expression (1 - sal ! 77). where sal is the salinity in ppt. Based on this result, 
differences between the open ocean salinity and the laboratory tank salinity should not 
cause the predicted spectral levels to be underestimated by more than about 0.1 dB. 

Temperature effects, however, have a significant impact on the sound 
energy radiated by large raindrops. Jacobus (1991) determined that the energy 
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radiated by large raindrops was linearly correlated with the absolute temperature 
difference between the raindrop and the water surface. For example, for 4.2 mm 
drops in fresh water Jacobus obtained the expression E - 540 + 38 -AT, where E is 
the raindrop acoustic energy in picojoules, and AT the drop-surface temperature 
difference in °C. Because Jacobus's spectral analysis technique was flawed, this 
expression is not strictly correct. However, the expression can be used to estimate the 
effect of drop-surface temperature differences if it is assumed his spectral analysis 
error affected all of his results uniformly. Using Jacobus's expression, raindrops of 
4.2 mm diameter whose temperature difference with the surface is 20° C will radiate 
approximately 3.8 dB more energy than raindrops that are at the same temperature as 
the surface. No measurements of raindrop temperatures have been made at the OTP, 
but it is possible that drop-surface temperature differences at sea were significantly 
larger than for the laborator)’ measurements. 

The effect of surface tension variations on the sound radiated by large 
drops has not yet been studied. The surface tension of the tank water was lower than 
for clean seawater by as much as 20 dynes/cm for some of the measurements, 
suggesting that surface contamination may have been present. A study of the surface 
tension effects on the sound generation by large raindrops should be performed. 
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Another factor that would tend to make the predicted spectral levels too 
low is that Type I raindrop bubbles were excluded from the prediction. As mentioned 
previously, the spectral peak at 15 kHz that is associated with the Type I raindrop 
bubbles is absent in the presence of moderate and heavy rainfall, indicating that Type 
I bubble formation is diminished. Some unknown fraction of small raindrop bubbles 
may still be present during moderate and heavy rainfall. 

A study currently in progress of the damping constants of raindrop 
bubbles in fresh and salt water has revealed that the lower frequency bubble signals 
collected in the redwood tank are contaminated to some extent by scattering of sound 
from the tank walls. It was suspected that tank scattering would increase the energy 
measured for individual raindrops, which in turn would result in higher than actual 
spectral level predictions. To estimate the magnitude of the error due to tank 
reverberation, dominant bubble signals of various frequencies were compared with the 
energies of simulated bubble signals with the same amplitude and frequency, but w'ith 
damping constants equal to the theoretical values. The energy difference between the 
measured and theoretical signals was on the order of ± 1 dB, with no observed 
dependence on frequency. 
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C. SUGGESTIONS FOR FUTURE WORK 



As mentioned earlier, a statistical approach could be used to obtain an inversion 
matrix A. This would be accomplished by making simultaneous measurements of 
raindrop size distributions and rainfall spectral levels. For each frequency bin of the 
rainfall spectrum, a multiple linear regression could then be performed between the 
rainfall spectral level and the drop rate densities for the various drop sizes. The 
regression coefficients obtained for a single frequency bin would correspond to one 
row of the A matrix described earlier. 

Comparing the A matrix coefficients obtained asing the statistical approach 
with the similar coefficients obtained from the laboratory measurements would shed 
some light on the physical processes of raindrop underwater sound generation. 
Variations in the coefficients due to wind speed and surface roughness could also be 
studied. The procurement of a raindrop size measuring device for data collection at 
the Ocean Test Platform site is anticipated. 
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VI. CONCLUSIONS 



The original goal of this study was to develop an inversion technique for 
measuring rainfall using laboratory measurements of the underwater sound generated 
by large raindrops. The laboratory measurements conducted here provide the basis 
for the for\\'ard problem required as part of the inversion. Differences between the 
measured and predicted rainfall spectra need to be accounted for before inversion of 
the rainfall spectra to obtain rainfall rate is attempted. The work accomplished here 
provided additional insights into the physical mechanisms for underwater sound 
generation by individual raindrops. For instance, the study revealed that; 

• The far field impact sound of rainfall consists primarily of impulsive broadband 
noise. 

• The energy of the large raindrop bubble sound tends to become smaller and 
more broadband in nature for decreasing drop size. 

• The dominant bubbles produced by the drop splash are responsible for most of 
the underwater sound energy of individual large raindrops. 

A statistical approach for relating underwater spectral levels to raindrop size 
distributions would be fruitful, and may reveal more about the physical processes of 
underwater sound generation by rainfall. 
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APPENDIX A 



DERIVATION OF THE EQUATION FOR THE ENERGY SPECTRUM 

The continuous Fourier Transform for aperiodic signals is given by: 

X{f) - 

Assume the signal is discretized using t - nAt, where n ranges from -NH to NH (in 
the limit as A^goes to infinity). The discrete form of the above equation is then: 

V2 -Jk2%nSt V2 - /2 n X /? 

X(f.) = Y. x{n)e~^^ At - AlY ^(n)e ^ - At • X{k) 

- V2 - V2 

where XiT^ ) becomes At ■ X{k ). and X{k ) is the Discrete Fourier Transform of the 
signal .V . 

By Rayleigh’s Theorem for an aperiodic signal, the energy is conserved per the 
follovring equation; 



E - JLy(/)IV/ - ^\X{f)\-df 



( 37 ) 



which, after substituting the discrete forms of x{t), X{f), and df - MNAt, becomes; 
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(38) 



-NI2 -K2 -N!2 NA I 

or, after canceling like terms; 

E ~ Y. - — E 

-A"2 N _,V2 

The units of ;r(n ) are Pa (the 1 m on-axis pressure), while the units oi X{k) are Pa 
per bin. Equation 39 can be rewritten so that the Fourier Spectrum is summed over all 
frequencies, where the frequency resolution is given by df - M(NAt ) : 

V2 

E = E \Mkf-df 
N-df^2 

The individual terms on the right hand side of Equation 40 represent the energy 
spectrum of the finite duration pressure signal, in units of Pa"- s/Hz. Since the 
frequency spectrum is symmetric, the negative frequencies can be mapped into the 
positive frequencies b>' multiplying the spectral levels by two, yielding: 

E{f) - • \X{k)\\ 0 (41) 

N-df 

where E{f) represents the energy density spectral level in Pa^- s/Hz for the frequency 
f - k • df as the index k varies from 0 to Nil. 
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To verify Equation 41, a damped sinusoid given by /?(0 - - 008 ( 211 /^ /) was 

used to simulate an actual bubble signal. The value of a was taken to be 1/(5 •10'^) s ', 
and w as assumed to be 5000 Hz. The signal was sampled at an interval of 1 ^lsec, 
with a total record length of 2048 points. A plot of the signal is shown in Figure 71. 
The energy density spectrum of the signal was then calculated using Equation 41 . The 
resulting spectrum is shown in Figure 72. The total signal energy was then calculated 
in both the time domain and the frequency domain to yield the following results: 

n - S - 1 

52 Pinfdt - 1.2446-10-^ Pa'-s 

n-O (42) 

52 E{f)df- 1.2495-10'" Pa^-s 

r-o 

The resulting values are in good agreement. 
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Figure 71. Simulated bubble signal for spectrum verification. 
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Figiire 72. Energy density specinim of simulated bubble signal. 
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APPENDIX B 



DERIVATION OF THE BUBBLE ENERGY EQUATION 

Let (/, r) be the instantaneous on-axis acoustic pressure of a dipole source at 
a range r, and p (/, r, 0 ) be the pressure at the range r and angle 6. For a dipole 
source, p{t, r, 6 ) - r) cos d, where 6 is the angle with the dipole axis. The 
radiated acoustical energy E can be written as: 

E - \ INTENSITY) dA dt - \ \{p-u)dAdt (43) 

1 A i A 

where u is the particle velocity and A is the surface area. The integration in area is 
performed over the hemisphere from the dipole axis to the surface. 

When k ■ R >3, the relationship u = p t {pc) can be used with less than 10% 
error in magnitude, where pc is the acoustic impedance and k is the wavenumber. 
Performing a surface integration with a ring element given by (In • r • sin 6 ) rdd 
and substituting the expression for p in terms of p^, we have: 
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^ 2 

E - ( ^ ^ ) r f d%\nd (id dl 

'ro 
,2 71 



,Zn r f 2 . 

(- ) dt 

3oc 



which can be expressed in digitized form as: 



r- T 



A/ 

3 p c f-o 



where T is the time duration of the signal (eg., the record length). 



(44) 



(45) 
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APPENDIX C 



DERIVATION OF THE RAINFALL SPECTRUM EQUATION 

Given DRD {d ) (the drop rate density for a given raindrop size in units of 
drops/m\s), and E{f.d) (the ] meter on-axis energy density spectral values for a given 
drop size in units of Pa^* s/Hz), the spectral values at a depth due to raindrops of a 

given size can be determined by summing the individual contributions from the bubble 
and impact dipole sources distributed over the surface. 

The spectral values at the depth must be obtained by integrating the 
contributions to the underw ater sound from each drop size over the surface area, while 
accounting for the dipole nature of the bubble and impact sound. The spectral values 
for each drop size can then be summed to yield the overall rainfall spectrum. The 
geometr>- of the problem is illustrated in Figure 73. 

The pressure /^at a range r and angle 6 from a surface dipole is equal to the on- 
axis pressure times the ratio cos(0 )*r^ /r, w'here is 1 meter. If the range r is such 
that the plane wave approximation is valid, then the intensity is approximately equal 
to /7* / {pc), where pc is the acoustic impedance. If the intensity is proportional to the 
square of the pressure, then the intensity at the depth z^ due to an arbitrary dipole will 
be equal to the on-axis intensity of the dipole at one meter times the square of the ratio 
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cos(d )T, tr. The efTecl of attenuation has been neglected here. Accounting for the 



dipole variation of intensity, the total rainfall spectrum at a depth can be expressed 
as: 



where the area integration is carried out over the ocean surface. Referring to Figure 
73, the surface area element dA for a ring integration is given hy dA - In • h • dh , 
where h varies from zero to r • sin iff (f). The angle f (f) corresponds to the 
beamwidth at a frequency T for an upward looking receiver at the depth z^. Equation 
46 can be transformed to an integral in d using the substitutions z^ = rcos 6, 
h = rsin 6. and dh = z, sec^ 6 dd . which leads to: 

o 




RS{f) = f E{f.d)- DRD{d)-co^d’^m6dd] (^”7) 



^ 0 



0 



Finally, evaluation of the integral leads to: 



RS{r) - Y,T^-^^^^^)-DRD{d)-s\n^iff{f)-\m^ 



(48) 



d 



where rf has been replaced by 1 m^. 
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Figure 73. Geomeirv' for ihe ring iniegralion of the surface sound field. 
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